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ABSTRACT 

Wc present a catalog ol 68 galaxy clusters, of which 19 are new discoveries, detected via the Sunyaev- 
Zel'dovich effect (SZ) at 148 GHz in the Atacama Cosmology Telescope (ACT) survey on the celestial 
equator. With this addition, the ACT collaboration has reported a total of 91 optically confirmed, 
SZ detected clusters. The 504 square degree survey region includes 270 square degrees of overlap 
with SDSS Stripe 82, permitting the confirmation of SZ cluster candidates in deep archival optical 
data. The subsample of 48 clusters within Stripe 82 is estimated to be 90% complete for M 50 o c > 
4.5 x 1O 14 M and redshifts 0.15 < z < 0.8. While a full suite of matched filters is used to detect the 
clusters, the sample is studied further through a "Profile Based Amplitude Analysis" using a statistic 
derived from a single filter at a fixed $500 = 5.' 9 angular scale. This new approach applies a one- 
parameter family of "Universal Pressure Profiles" (UPP) with associated scaling laws, derived from 
X-ray measurements of nearby clusters, to break the degeneracy between the cluster characteristic size 
(i?50o) and the integrated Compton parameter (Y 500 ). Three additional models of cluster physics are 
used to investigate a range of scaling relations beyond the UPP prescription. Assuming a concordance 
cosmology, the UPP scalings are found to be nearly identical to an adiabatic model, while a model 
incorporating non-thermal pressure better matches dynamical mass measurements and masses from 
the South Pole Telescope. A complete, high signal to noise ratio subsample of 15 ACT clusters 
is used to obtain cosmological constraints. We demonstrate, using fixed scaling relations, how the 
constraints depend on the assumed gas model if only SZ measurements are used, and show that 
constraints from SZ data are limited by uncertainty in the scaling relation parameters rather than 
sample size or measurement uncertainty. We next add in seven clusters from the ACT Southern 
survey, including their dynamical mass measurements, which are based on galaxy velocity dispersions 
and thus are independent of the gas physics. In combination with WMAP7 these data simultaneously 
constrain the scaling relation and cosmological parameters, yielding 68% confidence ranges described 
by erg = 0.829 ± 0.024 and TO = 0.292 ± 0.025. We consider these results in the context of constraints 
from CMB and other cluster studies. The constraints arise mainly due to the inclusion of the dynamical 
mass information and do not require strong priors on the SZ scaling relation parameters. The results 
include marginalization over a 15% bias in dynamical masses relative to the true halo mass. In an 
extension to ACDM that incorporates non-zero neutrino mass density, we combine our data with 
WMAP7, Baryon Acoustic Oscillation data, and Hubble constant measurements to constrain the sum 
of the neutrino mass species to be J2u m " < 0.29 °V (95% confidence limit). 

Subject headings: cosmology:cosmic microwave background - cosmology:obscrvations - galax- 
iesxlusters - Sunyaev-Zel'dovich Effect 
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1. INTRODUCTION 

Galaxy clusters are sensitive tracers of the growth of 
structure in the Universe. The measurement of their 
evolving abundance with redshift has the potential to 
provide constraints on cosmological parameters that are 
complementary to other measurements, such as the an- 
gular power spectrum of the cosmic microwave back- 
ground (CMB; e.g., Komatsu et al. 2011; Dunkley et al. 
2011; Keisler et al. 2011; Story et al. 2012), Type la su- 
pernovae (e.g., Hicken et al. 2009; Lampeitl et al. 2010; 
Suzuki et al. 2012), or baryon acoustic oscillations mea- 
sured in galaxy correlation functions (e.g., Percival et al. 
2010). 

There is a long history of using optical (e.g., Abell 
1958; Lumsden et al. 1992; Goto et al. 2002; Lopes et al. 
2004; Miller et al. 2005; Koester et al. 2007; Hao et al. 
2010; Szabo et al. 2011; Wen et al. 2009, 2012) and 
X-ray (e.g., Henry ct al. 1992; Bohringer et al. 2004; 
Burenin et al. 2007; Mehrtens et al. 2012) surveys to 
search for galaxy clusters. Data from such surveys 
offered an early indication of an £l m < 1 universe 
(e.g., Bahcall & Cen 1992). Recent results have demon- 
strated the power of modern optical and X-ray sur- 
veys for constraining cosmology (e.g., Vikhlinin et al. 
2009b; Mantz et al. 2010b; Rozo et al. 2010). A promis- 
ing method for both detecting clusters in optical sur- 
veys and simultaneously providing mass estimates is to 
use weak gravitational lensing shear selection, and the 
first such samples using this technique have recently ap- 
peared (e.g., Wittman ct al. 2006; Miyazaki et al. 2007). 
Within the last few years, cluster surveys exploiting 
the Sunyaev-Zel'dovich effect (SZ; Sunyaev & Zcl'dovich 
1970) have also begun to deliver cluster samples 
(e.g., Staniszcwski ct al. 2009; Marriage et al. 2011; 
Williamson et al. 2011; Planck Collaboration 2011a; 
Reichardt et al. 2012b) and constraints on cosmological 
parameters (Vanderlinde et al. 2010; Sehgal et al. 2011; 
Benson et al. 2011; Reichardt et al. 2012b). 

The thermal SZ effect is the inverse Compton scat- 
tering of CMB photons by electrons within the hot 
(~ 10 7 ~ 8 K) atmospheres of galaxy clusters. This leads 
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to a spectral distortion in the direction of clusters, with 
the size of the effect being proportional to the volume- 
integrated thermal pressure and thus, in the adiabatic 
scenario, the total thermal energy of the cluster gas. 
Accordingly, this is correlated with cluster mass (e.g., 
Bonamente et al. 2008; Marrone et al. 2012; Sifon et al. 
2012; Planck Collaboration 2012a). Since the SZ sig- 
nal is not diminished due to luminosity distance, it is 
nearly redshift independent; in principle SZ surveys can 
detect all clusters in the Universe above a mass limit 
set by the survey noise level (e.g., Birkinshaw 1999; 
Carlstrom et al. 2002). 

Although current SZ cluster samples are small 
in comparison to existing X-ray and optical cluster 
catalogs, they provide very powerful complementary 
probes because they are sensitive to the high mass, 
high redshift cluster population (e.g., Brodwin et al. 
2010; Foley et al. 2011; Planck Collaboration 2011b; 
Menanteau et al. 2012a; Stalder et al. 2012). Many 
studies (e.g., Hoyle et al. 2011; Mortonson et al. 2011; 
Hotchkiss 2011; Harrison & Coles 2012) have noted that 
the discovery of a sufficiently massive cluster at high red- 
shift would be a challenge to ACDM cosmology, and the 
approximately redshift independent, mass-limited nature 
of SZ surveys means that they are well suited to reveal 
such objects if they exist. 

In this paper we describe the results of a search for 
galaxy clusters using the SZ effect in maps of the celes- 
tial equator obtained by the Atacama Cosmology Tele- 
scope (ACT; Swetz et al. 2011). ACT is a 6m telescope 
located in northern Chile that observes the sky in three 
frequency bands (centered at 148, 218, and 277 GHz) si- 
multaneously with arcminute resolution. During 2008, 
ACT surveyed a 455 deg 2 patch of the Southern sky, 
centered on S = —55 deg, detecting a number of SZ clus- 
ter candidates of which 23 were optically confirmed as 
massive clusters (Menanteau et al. 2010; Marriage ct al. 

2011) . The Equatorial survey area, on which we report in 
this work, was chosen to overlap the deep (r w 23.5 mag) 
optical data from the Sloan Digital Sky Survey (SDSS; 
Abazajian ct al. 2009) Stripe 82 region (S82 hereafter; 
Annis et al. 2011). Optical confirmation of our SZ clus- 
ter candidates is reported in Menanteau et al. (2012b), 
using data from SDSS and additional targeted optical 
and IR observations obtained at Apache Point Observa- 
tory. All clusters have photometric redshifts, and most 
have spectroscopic redshifts from a combination of SDSS 
and new observations at Gemini South. X-ray fluxes 
from the ROSAT All Sky Survey confirm that this is a 
massive cluster sample. The overlap of the ACT survey 
with SDSS has also enabled stacking analyses which char- 
acterize the SZ-signal as a function of halo mass from op- 
tically selected samples (Hand ct al. 2011; Schgal ct al. 

2012) , as well as a first detection of the kinetic SZ effect 
from the correlation of positions and redshifts of lumi- 
nous red galaxies with temperature in the ACT maps 
(Hand et al. 2012). 

The structure of this paper is as follows. In Section 2 
we describe the processing of the ACT data used in this 
work and the cluster detection algorithm. In Section 3 
we describe our approach to relating the cluster signal 
in filtered SZ maps to cluster mass, and obtain mass 
estimates for our cluster sample. In Section 4 we com- 
pare our catalog to other SZ, optical, and X-ray selected 
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cluster catalogs. In Section 5 we obtain constraints on 
cosmological parameters using the ACT cluster sample. 
In an Appendix, we present analogous SZ signal mea- 
surements and mass estimates for ACT's Southern field 
clusters, using deeper data obtained over the course of 
the 2009-2010 observing seasons. 

Where it is necessary to adopt a fiducial cosmol- 
ogy, we assume £l m = 0.3, Q\ = 0.7, and Hq = 
7O/170 km s _1 Mpc" 1 (hro = 1), unless stated oth- 
erwise. Throughout this paper, cluster mass is mea- 
sured within a characteristic radius with respect to the 
critical density such that, e.g., M^qqc is defined as the 
mass measured within the radius (-R500) at which the 
enclosed mean density is 500 times the critical density 
at the cluster redshift. The function E{z) denotes the 
evolution of the Hubble parameter with redshift (i.e., 
E{z) = [O m (l + z) 3 + Oa] 1 / 2 for a universe with il k = 
and negligible radiation density). Uncertainties and er- 
ror bars are specified at the 1-a level. 

2. MAPS AND CLUSTER DETECTION 

In this section we discuss the detection of galaxy clus- 
ters in the ACT Equatorial maps at 148 GHz. The maps 
are filtered to enhance structures whose shape matches 
the Universal Pressure Profile of Arnaud et al. (2010). 
The final cluster catalog consists of SZ candidates that 
have been confirmed in optical or IR imaging. 

2.1. Equatorial Maps 

ACT's observations during the 2009 and 2010 seasons 
were concentrated on the celestial equator. For this study 
we make use of the 504 square degree deep, contiguous 
region spanning from 20 h 16 m 00 s to 3 h 52 m 24 s in right as- 
cension and from — 2°07' to 2°18' in declination. This 
region includes 270 square degrees of overlap with S82, 
which extends to 20 h 39 m in R.A. and ±1°15' in declina- 
tion. As shown in Figure 1, the S82 region corresponds 
to the lowest-noise region of the ACT Equatorial maps. 

The bolometer time-stream data are acquired while 
scanning the telescope in azimuth at fixed elevation. 
Cross-linked data are obtained by observing the same 
celestial region at two telescope pointings that produce 
approximately orthogonal scan directions. Because the 
scan strategy was optimized for simultaneous observa- 
tion with the 148 GHz and 218 GHz arrays (the centers 
of which are separated by approximately 33' when pro- 
jected onto the sky), the regions beyond declinations of 
±1°40' arc not well cross-linked. 

The ACT data reduction pipeline and map-making 
procedure are described in Diinner et al. (2013), in the 
context of ACT's 2008 data. The time-stream bolometer 
data are screened for pathologies and then combined to 
obtain a maximum likelihood estimate of the microwave 
sky map (with 0.5' pixels) for each observing season. As 
a result of the cross-linked scan strategy, and the careful 
treatment of noise during map making, the ACT maps 
are unbiased at angular scales £ > 300. 

Due to realignments of the primary and secondary mir- 
rors, the telescope beams vary slightly between seasons 
but are stable over the course of each season. The tele- 
scope beams are determined from observations of Saturn, 
using the method described in Hincks et al. (2010), but 
with additional corrections to account for « 6" RMS 



pointing variation between observations made on differ- 
ent nights (Hassclfield et al., in prep.). The effective 
beam for the 148 GHz array differs negligibly between 
the 2009 and 2010 seasons, with a FWHM of 1.4' and 
solid angle (including the effects of pointing variation) of 
224 ± 2 nsr. 

Calibrations of ACT observations are based on fre- 
quent detector load curves, and atmospheric opacity wa- 
ter vapor measurements (Diinner et al. 2013). Absolute 
calibration of the ACT maps is achieved by compar- 
ing the large angular scale (300 < £ < 1100) signal 
from the 2010 season maps to the WMAP 95 GHz 7- 
yearmaps (Jarosik et al. 2010). Using a cross-correlation 
technique as described in Hajian et al. (2011), an abso- 
lute calibration uncertainty of « 2% is achieved. The 
inter-calibration of 2009 and 2010 is measured through 
a similar cross-correlation technique, with less than 1% 
error. 

2.2. Gas Pressure Model 

At several stages in the detection and analysis we will 
require a template for the intracluster gas pressure pro- 
file. To this end we adopt the "Universal Pressure Pro- 
file" (UPP) of Arnaud et al. (2010, hereafter A10), which 
includes mass dependence in the profile shape and has 
been calibrated to X-ray observations of nearby clusters. 
In this section we review the form of the UPP, and ob- 
tain several approximations that will be used in clus- 
ter detection (Section 2.3) and cluster property recovery 
(Section 3). 

In A10, the cluster electron pressure as a func- 
tion of physical radius r is modeled with a gen- 
eralized Navarro-Frenk- White (GNFW) profile 
(Nagai, Kravtsov, & Vikhlinin 2007), 

P(x) = Po (csco^T 7 (1 + (c 5 oo^) Q ) (7 " /3)/Q , (1) 

where x = r/R^o and Pq, C500, 7, a, (3 are fit parameters. 
The overall pressure normalization, under assumptions 
of self-similarity (i.e., the case when gravity is the sole 
process responsible for setting cluster properties), varies 
with mass and redshift according to 

P500 = [1.65 x 10" 3 ^ keV cm" 3 ] m 2/3 E 8/3 (z), (2) 

where E(z) is the ratio of the Hubble constant at redshift 
z to its present value, and 

to = M 500c /(3 x 10 14 ft^Me) (3) 

is a convenient mass parameter. Some deviation from 
strict self-similarity may be encoded via an additional 
mass dependence in the shape of the profile, yielding a 
form 

P(r) = P 500 m a ^p(x). (4) 

In this framework, A10 use X-ray observations of local 
(z < 0.2) clusters to obtain best-fit GNFW parameters 

[Po,c 50 o,7,a,/3] = [8-403 h% 2 , 1.177, 0.3081, 1.0510, 
5.4905], and an additional radial dependence described 
reasonably well by a p (x) = 0.22/(1 + 8x 3 ). 

Because hydrostatic mass estimates are used by A10 to 
assess the relationship between cluster mass and the pres- 
sure profile, there may be systematic differences when 
one makes use of an alternative mass proxy, such as 
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weak lensing or galaxy velocity dispersion. Simulations 
suggest that hydrostatic masses are under-estimates of 
the true cluster mass (e.g., Nagai et al. 2007). How- 
ever, there is little consensus among recent studies which 
compare X-ray hydrostatic and weak lensing mass mea- 
surements. For example, Mahdavi et al. (2012) find 
that hydrostatic masses are lower than weak lensing 
masses by about 10% at i?soo; Zhang et al. (2010) and 
Vikhlinin et al. (2009a) find reasonable agreement; while 
Planck Collaboration (2012a) find hydrostatic masses to 
be about 20% larger than weak lensing masses. There- 
fore in our initial treatment of the UPP we neglect this 
bias; later we will address this issue by adding degrees of 
freedom to allow for changes in the normalization of the 
pressure profile. 

The thermal SZ signal is related to the optical depth for 
Compton scattering along a given line of sight. For our 
pressure profile, and in the absence of rclativistic effects, 
this Compton parameter at projected angle 9 from the 
cluster center is 



y(0) 



0~T 

m e c' 



J ds p (y s 2 + (r. 



soop/psoo,) 



(5) 



where #500 = -R500/ Da{z) with Da(z) the angular di- 
ameter distance to redshift z, ot is the Thomson cross 
section, m c is the electron mass, and the integral in s 
is along the line of sight. Relativistic effects change this 
picture somewhat, but for convenience we will use the 
above definition of y{9) and apply the rclativistic correc- 
tion only when calculating the SZ signal associated with 
a particular y. 

To simplify the expression for the cluster pressure pro- 
file, we first consider the mass parameter m = 1 and 
factor the expression in equation (5) to get 



y(9,m = 1) = 10 A ° E(z) 2 t(9/9 500 ) 



(6) 



where t(x) is a dimcnsionlcss profile normalized to 

r(0) = 1, and 10 A ° = 4.950 x lO" 5 /^ 2 gives the nor- 
malization. The deviations from self-similarity are weak 
enough that we may model the changes in the profile 
shape with mass as simple adjustments to the normal- 
ization and angular scale of the profile. For the masses 
of interest here (1 < m < 10) we obtain 

y(0,m) w 10 A "E(z) 2 m 1+Bo T(m Co 9/9 5Q0 ) (7) 

with Bq = 0.08 and C = -0.025. This approx- 
imation reproduces the inner signal shape extremely 
well, with deviations increasing to the 0.5% level by 
$500- For 0.16*500 < 9 < 39 5no , the enclosed signal 

(j Q e d9' 2ii9'y{9' , m)) differs by less than 1% from the re- 
sults of the full computation. This parametrization of the 
cluster signal in terms of a normalization and dimension- 
less profile is not used for cluster detection (Section 2.3), 
but will motivate the formulation of scaling relations and 
permit the estimation of cluster masses (Section 3.1). 

The observed signal due to the SZ effect is a change in 
radiation intensity, expressed in units of CMB tempera- 
ture: 



AT (6) 



fsz y{6). 



(8) 



TCMB 

In the non-relativistic limit, the factor /sz depends only 



on the observed radiation frequency. Integrating this 
non-relativistic SZ spectral response over the nominal 
148 GHz array band-pass, we obtain an effective fre- 
quency of 146.9 GHz (Swctz et al. 2011). At this fre- 
quency, the formulae of Itoh et al. (1998) provide a spec- 
tral factor, including rclativistic effects for gas tempera- 
ture T e , of / sz (t) = -0.992/ re i(t) where t = k B T e /m e c 2 
and / re i(i) = 1 + 3.79t - 28. 2t 2 . This results in a 6% 
correction for a cluster with T = lOkeV. We use the 
scaling relation of Arnaud et al. (2005), t = -0.00848 x 
(m£(z)) -0 ' 585 , to express the mean temperature depen- 
dence in terms of the cluster mass and redshift. This 
yields a final form, fsz(m, z) = — 0.992/ rG i(m, z), which 
we use in all subsequent modeling of the SZ signal. 

2.3. Galaxy Cluster Detection 

In addition to the temperature decrements due to 
galaxy clusters, the ACT maps at 148 GHz contain con- 
tributions from the CMB, radio point sources, dusty 
galaxies, and noise from atmospheric fluctuations and 
the detectors. To detect galaxy clusters in the ACT maps 
we make use of a set of matched filters, with signal tem- 
plates based on the UPP through the integrated profile 
template t(9/9 5 oo)- 

We consider signal templates Sg 500 (9) = t(9/9 50 o) for 
#500 = 1'18 to 27' in increments of 1 1 8 . Each fixed 
angular scale corresponds to a physical scale that varies 
with redshift, but can be computed for a given cosmology. 
For each signal template we form an associated matched 
filter in Fourier space 



*fl S oo(*0 



1 B(k)S 9500 (k) 



£0™ 



N(k) 



(9) 



where B(k) is the product of the telescope beam re- 
sponse with the map pixel window function, N(k) is the 
(anisotropic) noise power spectrum of the map, and T,g 500 
is a normalization factor chosen so that, when applied 
to a map containing a beam-convolved cluster signal 
— AT[Sg 500 * B]{9) (in temperature units), the matched 
filter returns the central decrement —AT. 

Since the total power from the galaxy cluster SZ sig- 
nal is low compared to the CMB, atmospheric noise, and 
white noise that contaminate the cluster signal, we esti- 
mate the noise spectrum N(k) from the map directly. 
Bright (signal to noise ratio greater than five) point 
sources are masked from the map, with the masking ra- 
dius ranging from 2' for the dimmest sources to 1° for the 
brightest source. A plane is fit to the map signal (weight- 
ing by the inverse number of samples in each pixel) and 
removed, and the map is apodized within 0.2° of the map 
edges. 

For ACT, the effect of the noise term in the matched 
filter is to strongly suppress signal below £ 3000 (corre- 
sponding to scales larger than 7'). When combined with 
the signal template (and beam), the filters form band- 
passes centered at I ranging from roughly 2500 to 5000. 

The azimuthally averaged real space filter kernel cor- 
responding to #500 = 5f 9 is shown in Figure 2, and com- 
pared to both the ACT 148 GHz beam and the signal 
template Sg 500 . 

The true noise spectrum may vary somewhat over the 
map due to variations in atmospheric and detector noise 
levels, and thus the matched filter #e 50 o might be said 
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Figure 1. The portion of the ACT Equatorial survey region considered in this work. It spans from 20 h 16 m 00 s to 3 h 52 m 24 s in R.A. and 
from — 2°07' to 2° 18' in declination for a total of 504 square degrees. The overlap with Stripe 82 (dashed line) extends only to 20 h 39 m in 
R.A. and covers ±1°15' in declination, for a total of 270 square degrees. Circles identify the optically confirmed SZ-sclectcd galaxy clusters, 
with radius proportional to the signal to noise ratio of the detection (which ranges from 4 to 13). The gray-scale gives the sensitivity (in 
CMB uK) to detection of galaxy clusters, after filtering, for the matched filter with 8500 = 5'.9 (see Section 2.3). Inside the Stripe 82 region 
the median noise level is 44 uK, with one quarter of pixels having noise less (respectively, more) than 41 uK (46 |jK). Outside Stripe 82, 
the median level is 54 \xK, with one quarter of pixels having less (more) than 47 (64) uK noise. The higher noise, A'-shaped regions are 
due to breaks in the scan for calibration operations. 
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Figure 2. The azimuthally averaged real space matched filter ker- 
nel, proportional to ^ 5 'g(^)> f° r signal template with S500 = 5f9. 
Shown for reference are the ACT 148 GHz beam, and the cluster 
signal template £g/ g (0). While filters tuned to many different an- 
gular scales are used for cluster detection (Section 2.3), the 5'.9 filter 
is used for cluster characterization and cosmology (Section 3.1). 



to be sub-optimal at any point. The filter remains un- 
biased, however, and a reasonable estimate of the signal 
to noise ratio (S/N) may still be obtained by recognizing 
that the local noise level will be highly correlated with 
the number of observations contributing to a given map 
pixel. 

Prior to matched filtering, the ACT maps are condi- 
tioned in the same way as for noise estimation, except 
that the point sources are subtracted from the maps in- 
stead of being masked. (Subsequent analysis disregards 
regions near those point sources, amounting to 1% of the 
map area.) With the application of each filter we obtain 
a map of AT values. A section of a filtered map (for 
#500 = 5. '9) is shown in Figure 3. 
We characterize the noise in each filtered map by 
8 10 modeling the variance at position x as <J 2 {x) — Oq + 

f7 hits/ n hits(^) ! where nhits(3>) is the number of detector 
samples falling in the pixel at x. We fit constants an and 
Chits by binning in small ranges of 1/rihits- The fit is it- 
erated after excluding regions near pixels that are strong 
outliers to the noise model. Typically such pixels are 
near eventual galaxy cluster candidates. Figure 1 shows 
the noise map for the 9 500 = 5'9 filter. 

After forming the signal to noise ratio map 
—AT(x)/a(x), cluster candidates arc identified as all 
pixels with values exceeding four. The catalog of cluster 
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Figure 3. Section of the 148 GHz map (covering 18.7 deg 2 ) match-filtered with a GNFW profile of scale 6500 = 5'9. Point sources are 
removed prior to filtering. Three optically confirmed clusters with S/N > 4.9 are highlighted (see Table 7). Within this area, there are an 
additional 11 candidates (4 < S/N < 4.9) which are not confirmed as clusters in the SDSS data (and thus may be spurious detections or 
high-redshift clusters) . 



candidates contains positions, central decrements (AT), 
and the local map noise level. Candidates seen at multi- 
ple filter scales are cross-identified if the detection posi- 
tions are within 1'; the cluster candidate positions that 
we list come from the map where the cluster was most 
significantly detected. We adopt the largest S/N value 
obtained over the range of filter scales as the detection 
significance for each candidate. 

For a given candidate, the S/N tends to vary only 
weakly with the filter scale. The reconstructed cen- 
tral decrement AT varies weakly above filter scales of 
#500 ~ 3', as may be seen for the most significantly de- 
tected clusters in Figure 4. As is discussed in Section 3.1, 
only the results from the 6500 = 5f 9 matched filter are 
used for inferring masses, scaling relations, and cosmo- 
logical results. The corresponding physical scale may be 
determined, as a function of rcdshift, based on cluster 
distance. 

2.4. Galaxy Cluster Confirmation 

The cluster candidates obtained from the 148 GHz map 
analysis are confirmed using optical and infrared imag- 
ing. A complete discussion of this process may be found 
in Menanteau et al. (2012b). For the purposes of this 
work, we briefly summarize the confirmation process and 
the redshift limits of the sample (which must be under- 
stood in order to derive cosmological constraints) . These 
limits differ according to the depth of the optical imaging 
available over a given part of the map. 

Most cluster candidates are confirmed through the 
analysis of SDSS imaging. The ACT Equatorial sur- 
vey is almost entirely covered by SDSS archival data 
(Abazajian et al. 2009), with a central strip designed to 
overlap with the deep optical data (r sa 23.5) in the S82 




Figure 4. Central decrement and signal to noise ratio as a func- 
tion of filter scale for the 20 clusters in S82 detected with peak 
S/N > 5. Top panel: Although the central decrement is a model- 
dependent quantity, the value tends to be stable for filter scales 
of 0500 > 3'. Bottom panel: On each curve, the circular point 
identifies the filter scale at which the peak S/N was observed. The 
vertical dashed line shows the angular scale chosen for cluster prop- 
erty and cosmology analysis. 6500 = 5/9. Despite the apparent gap 
near S/N m 6, the clusters shown represent a single population. 
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region (Annis et al. 2011), as shown in Figure 1. For each 
ACT cluster candidate with peak S/N > 4, SDSS images 
are studied using an iterative photometric analysis to 
identify a brightest cluster galaxy and an associated red 
sequence of member galaxies. For cluster candidates sat- 
isfying the confirmation criteria, redshifts are obtained 
from either a photometric analysis of the images, from 
SDSS spectroscopy of bright cluster members, or from 
targeted multi-object spectroscopic follow-up. The red- 
shift limit of cluster confirmation using SDSS data alone 
is estimated to be z « 0.8 within S82 and z « 0.5 outside 
of S82. High significance clusters that are not confirmed 
are investigated through targeted observations, under the 
assumption that they may be high rcdshift clusters. 

Within the S82 region, 49 candidates are confirmed, 
with 44 of these confirmations resulting from analysis 
of SDSS data only. Targeted follow-up of the high S/N 
candidates was pursued at the Apache Point Observa- 
tory, yielding five more confirmations, all at z > 0.9. 
The follow-up in the S82 region is deemed complete to a 
S/N of 5.1, in the sense that all SZ candidates with ACT 
S/N > 5.1 have been targeted. It is thus this sample, 
and this region, that are considered for the cosmological 
analysis (in addition to a subset of the Marriage et al. 
(2011) sample; see Section 5). The purity of the S/N 
> 5.1 candidate list is high; in fact all candidates were 
confirmed as clusters. The completeness within S82, as 
a function of redshift, is estimated in Section 3.6. 

Outside of S82, 19 clusters are confirmed using SDSS 
DR8 data. High significance SZ detections in this region 
that are not confirmed in the DR8 data constitute good 
candidates for high redshift galaxy clusters and are being 
investigated in a targeted follow-up campaign. 

In Table 7 we present the catalog of 68 confirmed clus- 
ters. For each object we fist its coordinates, rcdshift (sec 
Mcnanteau et al. 2012b, for details), S/N of the detec- 
tion (we adopt the maximum S/N across the range of 
filters used), and SZ properties. Figure 5 shows postage 
stamp images of some high-significance clusters, taken 
from the filtered ACT 148 GHz maps. 

3. RECOVERED CLUSTER PROPERTIES 

In this section we develop a relationship between clus- 
ter mass and the expected signal in the ACT filtered 
maps. The form of the scaling relationship between the 
SZ observable and the cluster mass is based on the UPP, 
and parameters of that relationship are studied using 
models of cluster physics and dynamical mass measure- 
ments. We obtain masses for the ACT Equatorial clus- 
ters assuming a representative set of parameters. 

3.1. Profile Based Amplitude Analysis 

Scaling relations between cluster M^qq c mass and clus- 
ter SZ signal strength are often expressed in terms of 
bulk integrated Compton quantities, such as Y" 50 o, which 
are expected to be correlated to mass with low intrinsic 
scatter (e.g., Motl et al. 2005; Reid & Spergel 2006). To 
measure this quantity from SZ signal maps requires an 
estimate of the physical cluster scale (e.g., R500 1 ), which 
is difficult to constrain from the 148 GHz data alone. 

To obtain Y500 measurements, the Planck 
team has used X-ray luminosity based masses 

1 Msooc = (4tt/3) x 500 Pc {z)R 3 500 ; 6» 50 o = R 5 oo/D A (z). 



(Planck Collaboration 2011a) as well as more detailed 
X-ray and weak lensing studies (Planck Collaboration 
2012a) to constrain i?5oo- 

In the absence of suitable X-ray or optical measure- 
ments of SZ clusters, authors have constructed scal- 
ing relations based on alternative SZ based statistics, 
such as the amplitude returned by some particular fil- 
ter (Sehgal et al. 2011), or the maximum S/N over some 
ensemble of filters (Vanderlinde et al. 2010). Recogniz- 
ing that the cluster scale is poorly constrained by the 
filter ensemble, recent work from the South Pole Tele- 
scope has included a marginalization over the results re- 
turned by the ensemble of filters (e.g., Story et al. 2011; 
Reichardt et al. 2012b). 

As an alternative approach, we adopt a single filter to 
use for all clusters, regardless of redshift, and then ac- 
count for the mismatch between the filter scale and the 
cluster scale in a self-consistent way. In particular, we 
choose the filter associated with 8 50 o = 5'.9, and denote 
it as vE'(fc). One benefit of this approach is that, by using 
the same filter, the statistic can be computed using any 
map, facilitating comparison between surveys. In order 
to interpret the statistic returned by this filter, we model 
clusters as being well described, up to some overall ad- 
justments to the normalization and mass dependence, by 
the UPP (see Section 2.2). 

An estimate of the cluster central Compton parameter, 
based only on the non-relativistic SZ treatment, is given 

by 

AT , 

yo = - / sz 1 (m = 0,z = 0), (10) 

Jcmb 

where fsz( m = 0,z = 0) = —0.992 as explained in Sec- 
tion 2.2. This "uncorrected" central Compton parameter 
is used in place of AT to develop an interpretation of the 
SZ signal. This quantity is uncorrected in the sense that 
it is associated with the fixed angular scale filter and docs 
not include a rclativistic correction. 

For a cluster with SZ signal described by equation (7), 
the value of yo that we would expect to observe by ap- 
plying the filter to the beam-convolved map is 

y = W Ao E(z) 2 m 1+B "Q{9 500 /m c ")f rol (m, z) (11) 
where 

Q(0) = J (02*( fc ) B ( fe ) / d 2 9' e id '- k T(9'/6). (12) 

is the spatial convolution of the filter, the beam, and the 
cluster's unit-normalized integrated pressure profile. We 
note that in this formalism, #500 = R^oo/Da^z) is de- 
termined by the cluster mass and the cosmology (rather 
than being some independent parameter describing the 
angular scale of the pressure profile). 

The response function Q{9) for the Equatorial clusters 
is shown in Figure 6. It encapsulates the bias incurred 
in the central decrement estimate due to a mismatch be- 
tween the true cluster size and the size encoded in the 
filter, for the family of clusters described by the UPP. 
While this bias is in some cases substantial {Q sw 0.3), 
the function Q(8) is not strongly sensitive to the details 
of the assumed pressure profile (as demonstrated in Sec- 
tion 3.3), and the assumptions underlying this approach 
are not a significant departure from other analyses that 
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Figure 5. Postage stamp images (30' on a side) for the 10 highest S/N detections in the catalog (see Table 7), taken from the filtered 
ACT maps. The clusters are ordered by detection S/N, from top left to bottom right, and each postage stamp shown is filtered at the 
scale which optimizes the detection S/N. Note that J2327.4— 0204 is at the edge of the map. The greyscale is linear and runs from -350 /iK 
(black) to +100 AtK (white). 
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Figure 6. The response function used to reconstruct the cluster 
central decrement as a function of cluster angular size (solid line) . 
At 05oo = 5'9, the filter is perfectly matched and Q = 1. At 
scales slightly above 5'. 9, Q > 1 because such profiles have high 
in-band signal despite being an imperfect match, overall, to the 
template profile. For the definition of Q, see Section 3.1. Dotted 
line shows analogous function computed under the assumption that 
the cluster signal is described by the Planck Pressure Profile (see 
Section 3.3). 

rely on a family of cluster templates to extract a cluster 
observable. 

Equation (11) thus relates yo to cluster mass and red- 
shift while accounting for the impact of the filter on clus- 
ters whose angular size is determined by their mass and 
redshift. This relationship can be seen in Figure 7. 

The essence of our approach, then, will be to filter the 
maps with ^>(k) and for each confirmed cluster obtain 
AT and its error. This is equivalent to measuring yo, 
which can then be compared to the right hand side of 
equation (11). If the cluster redshift is also known, then 
for a given cosmology the only free parameter in the ex- 
pression for j/o is the mass parameter, m. 2 

2 With yo and z measurements in hand, one could certainly 



We refer to this approach, where a family of pressure 
profiles is used to model the amplitude of a source in 
a filtered map, as "Profile Based Amplitude Analysis" 
(PBAA). The reason we expect this approach to be effec- 
tive in SZ cluster analyses (and especially in cosmological 
surveys) is that the effects of angular diameter distance, 
telescope beam, and the spatial filtering are modeled in a 
way that accounts for the (mass and redshift dependent) 
cluster angular scale. For a given cosmology, and having 
computed Q{&) based on the UPP, the parameters {Aq, 
Bo, and Co) of the scaling relation between yo and mass 
have a physical significance and can be verified through 
measurements of yo, redshift, and mass for a suitable set 
of clusters. 

In order to compare the predictions of the UPP based 
formalism to models and other data sets, we introduce 
a more general relationship relating cluster mass to the 
uncorrected central Compton parameter. We allow for 
variations in the normalization, mass dependence, and 
scale evolution through parameters A, B, and C and 
model yo as 



~ = iq a 0+ a E (zY(M/M pivot ) 



1 + B +B, 



Q 



1 + Z 

1.5 



c 



7500 



/m 



Co 



/ re l(m, z). (13) 



To abbreviate the argument to Q{&), we will often simply 
write Q(m,z). The exponents (Aq, Bo, Co) remain fixed 
to the UPP model values of equations (6) and (7), except 
where otherwise noted. For a given data set or model, 
-Mpivot will be chosen to reduce covariance in the fit val- 
ues of A and B. In Table 1 we present the fit parameters 
for various models and data sets discussed in subsequent 
sections. In order to compare fits from data sets with 

proceed to solve equation (11) to obtain a mass for each cluster. 
Because we are treating mass as one of the independent variables, 
however, such an approach would produce biased mass estimates; 
sec Section 3.2. 
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Figure 7. Prediction, based on the UPP, for cluster signal in a 
map match-filtered with #500 = 5f9, in units of uncorrected central 
Compton parameter yo and apparent temperature decrement — AT 
at 148 GHz (Section 3.1). Solid lines trace constant masses of, 
from top to bottom, M 500c = 10 15 ,7 X 10 14 ,4 X 10 14 , and 2 X 
10 14 h~gM Q . Dotted lines are for the same masses, but with the 
scaling relation parameter C = 0.5 to show the redshift sensitivity 
to this parameter. Above 2 « 0.5, the scaling behavior of the 
observable yo with redshift is stronger for higher masses because 
their angular size is a better match to the cluster template and 
the redshift dependence in Q does not attenuate the scaling of the 
central decrement, yo oc E(z) 2 , as much as it does for lower masses. 
The dashed lines correspond to S/N > 4 and S/N > 5.1, based on 
the median noise level in the S82 region. 

different A/ p i vot , we also compute the normalization ex- 
ponent A m associated with M p i vot = 3 x 10 14 /i^Mq for 
each data set. In these terms, the UPP model described 
by equation (11) corresponds to (A m , B,C) = (0,0,0). 

One should not expect to be able to compare yo values 
measured by different instruments, because the bias fac- 
tor Q depends upon the beam and noise profiles of each 
telescope. For clusters sufficiently similar to the UPP, 
however, one should be able to compare the "corrected" 
central Compton parameter 

yo = yo/ [frei(m, z)Q(m, z)] , (14) 

or any other quantity that can be derived directly from 
the UPP prescription. Indeed, one of the goals of this 
analysis is to present cluster properties that can be 
straightforwardly included in a different analysis. 

The uncorrected central Compton parameter measure- 
ments (yo) for the ACT Equatorial clusters are presented 
in Table 7. They will be used in subsequent sections to 
estimate cluster properties (such as corrected SZ quan- 
tities and mass) and to constrain cosmological parame- 
ters. For the Southern cluster sample, analogous mea- 
surements are presented in the Appendix. Between the 
Equatorial and Southern cluster samples, the ACT col- 
laboration has reported a total of 91 optically confirmed, 
SZ detected clusters. 

3.2. Cluster Mass and SZ Quantity Estimates 

Given measurements of cluster yo and redshift, one 
cannot naively invert Equations (11) or (13) to obtain 
a mass estimate. Because of intrinsic scatter, measure- 
ment noise, and the non-trivial (very steep) cluster mass 
function, the mean mass at fixed SZ signal yo will be 
lower than the mass whose mean predicted SZ signal is 
yo. The bias due to noise is often referred to as "flux 



boosting" and can be corrected in a Bayesian analysis 
that accounts for the underlying distribution of flux den- 
sities (Coppin ct al. 2005). The bias due to intrinisic 
scatter, however, is not restricted to the low significance 
measurements. Considering the population of clusters 
(at fixed redshift) in the (logm,logyo) plane, the locus 
(logm| log yo) (i.e., the expectation value of the log of 
the mass for a given central Compton parameter) lies 
at lower mass than (log yo | log to) . This phenomenon 
has been discussed in the context of galaxy cluster sur- 
veys by, e.g., Mantz et al. (2010a, see also the review by 
Allen, Evrard, & Mantz 2011). 

The mass of a cluster, however, can be estimated if 
one has an expression for the cluster mass function. We 
adapt the Bayesian framework of Mantz ct al. (2010a) 
to this purpose. The posterior probability of the mass 
parameter to given the observation yg b is 

P(m\y° h ) ex P(y° h \m)P{m) 

= (J dj% P®?\ffl)P(W\rn)\ P(m) (15) 

where yjj r represents the "true" SZ signal in the absence 
of noise, P(yo b |yo r ) is the distribution of yg b given y^ 
and the observed noise <5yg b , and P(to) is proportional to 
the distribution of cluster masses at the cluster redshift. 
The distribution P(y^\m) of the noise- free cluster signal 
yo is assumed to be log-normal about the mean relation 
given by Equation (13), i.e., 

log5g r ~JV(logifo(m,«);<7? t ) (16) 

with dint denoting the intrinsic scatter. 

We use the results of Tinker et al. (2008) to compute 
the cluster mass function, assuming the fiducial ACDM 
cosmology, with eg = 0.8. Scaling the mass function 
by the comoving volume element at fixed solid angle, we 
obtain dN(< m,z)/dz, the number of clusters, per unit 
solid angle and per unit redshift, that have mass less than 
to. The probability of a cluster in this light cone having 
mass to and redshift z may then be taken as P(to, z) cx 
d 2 N(< to, z) I dz dm. We account for redshift uncertainty 
by marginalizing the cluster mass function P(to, z) over 
the cluster's redshift error to obtain an effective P(m) at 
the observed cluster redshift. 

The marginalized masses obtained using Equation (15) 
are presented in Table 8. For the ACT Southern cluster 
sample, these masses are presented in the Appendix. In 
each case, masses are presented for the UPP scaling re- 
lation as well as for scaling relation parameters fit to SZ 
signal models (see Section 3.4) or dynamical mass data 
(see Section 3.5). 

A similar approach may be taken to estimate the true 
values of SZ quantities, given the observed quantities. In 
this case we are effectively only undoing the noise bias, 
while intrinsic scatter affects the underlying population 
function. We are interested in 

- P(ffi b \W) / dm P{yS\m)P{m). (17) 

The resulting probability distribution is used to obtain 
marginalized estimates of yo, 6*500 (which should be inter- 
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Figure 8. Example probability distributions for cluster mass (up- 
per panel) and SZ signal strength parametrized as a mass according 
to equation (13) (lower panel). The solid line PDF is the result of a 
direct inversion of the scaling relation described by equation (13). 
The corrected PDF (dashed line) is obtained by accounting for the 
underlying population distribution (bold line; arbitrary normaliza- 
tion). The correction is computed according to equation (15) for 
the upper panel, and according to equation (17) for the lower panel. 
Curves shown correspond to ACT-CL J0022.2-0036. 

prctcd as giving the scale of the pressure profile rather 
than the scale of the mass density profile), I500 (esti- 
mated within the SZ-inferred 6500) and Q. These quanti- 
ties are presented in Table 8, for the UPP scaling relation 
parameters. 

In Figure 8 we demonstrate the impact of the steep 
mass function on the inferred mass and SZ quantities. As 
the measurement noise decreases, the j/o measurements 
are less biased; but any intrinsic scatter in the yo~M 
relation will lead to bias in the naively estimated mass. 

3.3. The Planck Pressure Profile 

In Planck Collaboration (2012b), data for 62 mas- 
sive clusters from the Planck all-sky Early Sunyaev- 
Zcl'dovich cluster sample (Planck Collaboration 2011a) 
are analyzed to obtain a " Planck pressure profile" (PPP) 
based on measurements of the SZ signal. Integrating 
the PPP along lines of sight, the central pressure is 20% 
lower than the UPP but is higher than the UPP outside 
of 0.5i?5oo- Planck finds overall consistency between re- 
sults obtained with the UPP and the PPP. 

We assess the difference in inferred mass due to this 
alternative pressure profile by re-analyzing the ACT y$ 
using the PPP. A bias function Q is computed as in Equa- 
tion (12), but with the normalized Compton profile r as- 
sociated to the PPP (see Figure 6). We additionally 
determine scaling parameters, compatible with Equa- 
tion (7), of (10 Ao ,B ,C ) = (4.153 x 10~ 5 ^ 2 , 0.12, 0). 
Note that while the bias function shows increased sen- 
sitivity, compared to the UPP, for 8 50 o < 9', this is 
compensated for by the lower normalization factor 1Q A °. 
For the Equatorial cluster sample, we find the PPP 
masses to be well described by a simple mean shift of 
M.^P = 1.015 Mgj£, with 3% RMS scatter. Note that 
this is only a statement about the dependence of the ACT 
results on the assumed pressure profile; experiments that 
probe different angular scales may be more or less sensi- 



tive to such a change. 

While the change in inferred masses is in this case neg- 
ligible, we reiterate that our fully parametrized relation- 
ship between SZ signal and mass (Equation (13)) allows 
for freedom in the normalization, mass dependence, and 
evolution of cluster concentration with redshift. Mass or 
cosmological parameter estimation can be computed af- 
ter fixing these parameters based on any chosen pressure 
profile, model, simulation, or data set; all that is required 
is to compensate for the mismatch of our assumed pres- 
sure profile to the true mean pressure profile. 

3.4. Scaling Relation Calibration from SZ Models 

The previous sections have described a general ap- 
proach that relates cluster mass and redshift to SZ signal 
in a filtered map, given values for the scaling relation pa- 
rameters. In this section we obtain scaling relation pa- 
rameters based on three models for cluster gas physics. 
While the ACT data will be interpreted using each of 
these results, we do not yet consider any ACT data ex- 
plicitly. 

Current models for the SZ signal from clusters include 
contributions from non-thermal pressure support, star 
formation, and energy feedback and are calibrated to 
match detailed hydrodynamical studies and X-ray or op- 
tical observations (Shaw et al. 2010; Bode et al. 2012). 
Such models provide a useful testing ground for the as- 
sumptions and methodology of our approach to predict- 
ing SZ signal based on cluster mass. While models may 
suffer from incomplete modeling of relevant physical ef- 
fects, they are less vulnerable to some measurement bi- 
ases (e.g., by providing a cluster mass and alleviating the 
need for secondary mass proxies). In order to explore 
the current uncertainty in the SZ-mass scaling relation, 
we consider simulated sky maps based on three models 
of cluster SZ signal that include different treatments of 
cluster physics. 

Our study will center on maps of SZ signal produced 
from the SZ models and structure formation simula- 
tions of Bode et al. (2012, hereafter B12). The N- 
body simulations (Bode & Ostriker 2003) are obtained 
in a Trce-Particle-Mcsh framework, in which dark mat- 
ter halos have been identified by a friends-of-friends algo- 
rithm. The intracluster medium (ICM) of massive halos 
is subsequently added, following a hydrostatic equilib- 
rium prescription, and calibrated to X-ray and optical 
data (Bode, Ostriker, & Vikhlinin 2009). The density 
and temperature of the ICM of lower mass halos and 
the IGM are modeled as a virializcd ideal gas with den- 
sity (assuming cosmic baryon fraction fi;,/f2 m = 0.167) 
and kinematics that follow the dark matter. 

To complement the model of B12, we also consider 
the Adiabatic and Nonthermal20 models described in 
Trac et al. (2011), which make use of the same N-body 
results as B12. In the Adiabatic model the absence of 
feedback and star formation leads to a higher gas fraction 
than in the B12 model. In the Nonthermal20 model, 20% 
of the hydrostatic pressure is assumed to be nonthermal, 
leading to substantially less SZ signal compared to the 
B12 model. The SZ-mass relations derived from these 
two models are thus interpreted as, respectively, upper 
and lower bounds on the SZ signal. 

The model of B12 differs from those in Trac et al. 
(2011) through a more detailed handling of non-thermal 
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pressure support, which is tied to the dynamical state 
of the cluster and is allowed to vary over the clus- 
ter extent. Both B12 and the similar treatment of 
Shaw et al. (2010) make use of hydrodynamic simula- 
tions (Nagai et al. 2007; Lau et al. 2009; Battaglia et al. 
2012) to understand these non-thermal contributions. 

To calibrate our scaling relation approach to these 
models, we make use of light-cone integrated maps of 
the thermal and kinetic SZ at 145 GHz (constructed as 
in Sehgal et al. 2010), and the associated catalog of clus- 
ter positions and masses. A set of 192 non-overlapping 
patches of area 18.2 deg 2 each are extracted from the 
simulated map, and convolved with the ACT 148 GHz 
beam to simulate observation with the telescope. The 
maps arc then filtered with the same filter ty(k) that was 
used for the ACT Equatorial clusters. Because the filter- 
ing is a linear operation, it is counter-productive to the 
purpose of calibration and intrinsic scatter estimation to 
add noise (CMB, detector noise) to the simulated signal 
map, and so we do not. The uncorrected central decre- 
ments are extracted and used to constrain the parameters 
of Equation (13). To probe the high-mass regime, only 
the 257 clusters having M 500c > 4.3 x 10 14 h^M Q and 
0.2 < z < 1.4 are considered; the fit is performed around 
M p i vot = 5.5 x 10 14 /i7 O 1 M0. The intrinsic scatter of the 
relation is also obtained from the RMS of the residuals. 
For the B12 model, the residuals of the fit are plotted 
against mass in Figure 9. 

For each of the three models, fit parameters are pre- 
sented in Table 1. The mass dependence is consistent, 
in all cases, with the UPP prediction (B w 0), and ad- 
ditional redshift dependence is only present in the Non- 
thermal20 model. Only the Adiabatic model is consistent 
in its normalization with the UPP value. This is despite 
the explicit calibration, in B12, of the mean pressure pro- 
file to the UPP at i?soo- The origin of this inconsistency 
is due to the relative shallowness of the mean pressure 
profile in B12 compared to the UPP. Thus, the profiles 
in B12 have less total signal within -R500, where ACT is 
sensitive. 

The scaling relation parameters obtained for the Adi- 
abatic model are sufficiently close to zero (i.e., to the 
UPP scaling prediction), that we drop them from fur- 
ther consideration. While the B12 normalization lies 
somewhat below the UPP prediction, we note that even 
lower normalizations (such as that found in the Nonthcr- 
mal20 model) are favored by recent measurements of the 
SZ contribution to the CMB angular power spectrum 
(Dunkley et al. 2011; Reichardt et al. 2012a). 

We thus proceed to consider quantities derived from 
each of the UPP, B12, and Nonthcrmal20 scaling rela- 
tion parameter sets. Mass estimates for the B12 and 
Nonthcrmal20 models are computed for the ACT Equa- 
torial cluster sample as described in Section 3.2, and are 
presented alongside the UPP estimates in Table 8. 

3.5. Scaling Relation Calibration from Dynamical 
Masses 

Sifon et al. (2012, hereafter S12) measure galaxy ve- 
locity dispersions to obtain mass estimates for clusters in 
ACT's Southern field. S12 also present the uncorrected 
central Compton parameter measurements yo and the 
corrected versions yo obtained as described in Section 3.2 
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Figure 9. Residuals of the scaling relation fit for the B12 model 
(Section 3.4). Only clusters with 0.2 < z < 1.4 and M 5 oo c > 4.3 X 
10 14 (Iju'Mq (indicated by dotted line) are used for the fit. The 
scatter in the relation is measured from the RMS of the residuals. 



and presented in the Appendix. S12 perform power law 
fits of both the uncorrected (yo) and corrected (yo) cen- 
tral Compton parameters to the dynamical masses to 
establish scaling relations for those cluster observables. 

Here, we fit the full scaling relation of equation (13) 
to the dynamical mass data for the 16 z > 0.3 clusters 
from the Southern field that were detected by ACT and 
observed by S12. We do not use the scaling relation 
as estimated in S12, because the parametrization of the 
scaling relation in that study is different. Also, the lin- 
ear regression that is used in S12 (the bisector algorithm 
of Akritas & Bershady (1996)) is not suited to predict- 
ing the SZ signal given only the mass, which is the aim 
in formulating the cluster abundance likelihood in Sec- 
tion 5. 3 

Dynamical masses are estimated in SI 2 for each cluster 
based on an average of 60 member galaxy spectroscopic 
rcdshifts. For each cluster, the galaxy velocity dispersion 
Sbi is interpreted according to the simulation based re- 
sults of Evrard et al. (2008), who find that the dark mat- 
ter velocity dispersion odm is related to the halo mass 
M 2 ooc by 

( 0.7 xE(z) M 2QQc \ a 

where <r 15 = 1082.9 ± 4 km s" 1 and a = 0.3361 ± 0.0026. 
By inverting equation (18) and assuming that Ssi = 
(Jdm j S12 obtain dynamical estimates, which we will de- 
note by M^qqc, of the halo mass. As discussed in S12 and 
Evrard et al. (2008), the systematic bias between galaxy 
and dark matter velocity dispersions, b v = Sbi/cdm, is 
believed to be within 5% of unity. To account for this, 
and any other potential systematic biases in the dynam- 

3 The cluster abundance likelihood assumes a scaling relation 
where yo is the dependent variable and takes full account of the 
mass function; in this section we will use the likelihood-based ap- 
proach of Kelly (2007) which includes iterative estimation of the 
distribution of the independent variable. 
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Table 1 

Scaling relation parameters, fit to: (i) various SZ models (sec Section 3.4); (ii) the dynamical mass data of Sifon ct al. (2012) 
(Section 3.5); (iii) a cosmological MCMC including Southern and Equatorial cluster samples and dynamical mass data 
(Section 5.3). Scaling relation parameters A, B, and C arc defined as in equation (13), with Af p ; vot chosen to yield uncorrclatcd 

A and B. A m is the normalization parameter corresponding to Af p i vot = 3 X 10 14 h^ Mg and may be compared among rows. 

Parameters A m ,B and C indicate the level of deviation from the predictions based on the Universal Pressure Profile of 
Arnaud et al. (2010) (equations (6) and (7); shown for reference). The intrinsic scatter <Tj nt is defined as the square root of the 
variance of the observed logyo> in the absence of noise, relative to the mean relation defined by equation 13. The parameter C is 

fixed to when fitting scaling relations to dynamical masses. 



Description 




A 


Am 


B 


C 


fint 


Universal Pressure Profile (UPP) 















0.20 a 


Models (§3.4) 














B12 


5.5 


0.111 ±0.021 


-0.17 ±0.06 


-0.00 ±0.20 


-0.04 ±0.37 


0.20 


Nonthermal20 


5.5 


-0.003 ± 0.020 


-0.29 ±0.06 


0.00 ±0.20 


0.67 ±0.47 


0.21 


Adiabatic 


5.5 


0.241 ±0.020 


-0.02 ± 0.06 


-0.08 ±0.20 


0.10 ±0.43 


0.21 


Dynamical mass data (§3.5) 














All clusters 


7.5 


0.237 ±0.060 


-0.21 ±0.21 


0.03 ±0.51 





0.31 ±0.13 


Excluding J0102 


7.5 


0.205 ± 0.045 


-0.11 ±0.15 


-0.28 ±0.35 





0.19 ±0.10 


Full cosmological MCMC (§5.3) 














ACDM model 


7.0 


0.079 ±0.135 


-0.45 ±0.19 


0.36 ±0.36 


0.43 ±0.62 


0.42 ±0.19 


u>CDM model 


7.0 


0.065 ±0.153 


-0.46 ±0.21 


0.36 ±0.35 


0.34 ± 0.65 


0.45 ±0.20 



a This value, based on the B12 model value, is used for results computed for the UPP scaling relation parameters that also require 
a value for the intrinsic scatter. 



ical mass estimates, we introduce the parameter 

^dyn / M$&\ (1Q) 

P \M 200c / V ' 

Based on a velocity dispersion bias of b v = 1.00 ± 0.05, 
the equivalent mass bias is (3 dyn = 1.00 ± 0.15. For the 
present discussion, we disregard this bias in order to dis- 
tinguish its effects from other calibration issues. How- 
ever, in the cosmological parameter analysis of Section 5 
wc include /3 d y n as a nuisance parameter and discuss its 
impact on the cosmological parameter constraints. 

The 7/0 measurements associated with the Southern 
sample of clusters are obtained using a filter matched 
to the noise power spectrum of the Southern field maps 
used in Sifon et al. (2012). Thus, while the signal tem- 
plate is the same, the full form of the filter VP and the 
associated response function Q differ slightly from the 
ones used on the Equatorial data. We apply the same 
correction for selection bias that was used by SI 2, and 
denote the corrected values as y§ orr . 

To convert the dynamical masses to M^qqc values, we 
model the cluster halo with a Navarro-Frenk- White pro- 
file (Navarro et al. 1995) with concentration parameters 
and uncertainties obtained from the fits of Duffy et al. 
(2008). For the fit we use a pivot mass of 7.5 x 
10 14 JijqMq, and fix the parameter C to (otherwise 
the fit is poorly constrained). We use the likelihood- 
based approach of Kelly (2007) to fit for the intrin- 
sic scatter along with the parameters A and B, given 
measurement errors on both independent and dependent 
variables. The scatter is modeled, as before, as an addi- 
tional Gaussian random contribution to log yo relative to 
the mean relation (?/o| m , z )- 

The fit parameters are presented in Table 1. A sub- 
stantial contribution to the scatter in the dynamical mass 
fits comes from the exceptional, merging cluster ACT- 
CL J0102-4915 ("El Gordo," Menanteau et al. 2012a): 
when this cluster is excluded from the fit, the scatter 
drops to 0.19±0.10, which is more consistent with fits 
based solely on models. 



In Figure 10 we plot the cluster SZ measurements 
against the dynamical masses, along with the best fit 
scaling relation. The scaling relations from the UPP and 
from the parameters fit to the B12 and Nonthermal20 
models are also shown. While the fit parameters for the 
dynamical mass data are consistent with either the B12 
or Nonthcrmal20 models, the dynamical mass data lie 
well below the mean scaling relation predicted by the 
UPP. These results reinforce the need to consider a broad 
range of possible scaling relation parameters, within our 
framework based on the UPP. We note, however, that 
the possibility of a systematic difference between dynam- 
ical masses and other mass proxies must be considered 
when comparing the parameters obtained in this section 
to other results. 

3.6. Completeness Estimate 

In this section we estimate the mass, as a function 
of redshift, above which the ACT cluster sample within 
SDSS Stripe 82 (which we will refer to as the S82 sample) 
is 90% complete. We consider the S82 sample as a whole 
(S/N > 4), and also consider the subsample that has 
complete high redshift follow-up (S/N > 5.1). 

To obtain a model for the SZ signal in the maps as a 
function of mass for each filter scale, we apply the ap- 
proach of Section 3.1 to each matched filter ^e 500 i an d 
obtain a relationship between the uncorrected central 
Compton parameter yg 5aa and the cluster mass, as in 
Equations (11) and (13). Note that for the UPP, our 
formulation always has scaling parameters (A,B,C) ~ 
(0,0,0), but a matched bias function Qs 5OO (0) must be 
computed for each filter ^e 500 - 

For a cluster of a given mass and redshift, and for a 
particular filter scale, we estimate the average cluster de- 
tection probability by marginalizing over both the intrin- 
sic scatter in ye 50a as well as over the full space of noise 
realizations (taken from the noise map in the S82 re- 
gion). The largest detection probability, considering all 
filter scales, is taken as the detection probability. The 
resulting efficiency of detection as a function of mass 
and redshift is multiplied by the cluster mass function 
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Figure 10. Corrected central Compton parameter vs. dynami- 
cal mass for the 16 ACT-detected clusters presented in Sifon et al. 
(2012) for the Southern sample. Values on y-axis include factor of 
E(z)~ 2 , which arises in the derivation of yo in self-similar mod- 
els. The high signal outlier is "El Gordo" (ACT-CL J0102-4915, 
Menanteau et al. 2012a), an exceptional, merging system. The 
solid line represents the best fit of equation (13) with M p i vot = 
7.5 X 10 14 h~gM Q . The dashed line is for the fit with J0102-4915 
excluded. Dotted lines, from top to bottom, are computed for 
scaling relation parameters corresponding to the UPP, B12 and 
Nonthcrmal20 (z = 0.5) models. 

of Tinker et al. (2008) to calculate the mass above which 
each sample is 90% complete. 

The completeness mass levels are shown in Figure 11. 
Note that we also show results obtained for the B12 
scaling relation parameters. In this case we also ob- 
tained completeness estimates based in part on the fil- 
tered simulated maps. The central Compton parameters 
were measured in the filtered maps, and the S82 noise 
model was applied to generate a detection probability for 
each simulated cluster. Because of the small number of 
sufficiently high mass clusters in the model simulations, 
we have compensated for sample variance by reweight- 
ing the contribution of each cluster to correspond to the 
Tinker et al. (2008) mass function. 

In summary, the S82 sample is estimated to be 90% 
complete above masses of M 500c « 4.5 x 10 14 h^Mo, 
for 0.2 < z < 0.8. The S82 sample having S/N > 5.1 is 
estimated to be 90% complete, in 0.2 < z < 1.5, above 
M 50 oc ~ 5.1 x 10 14 H^qMq. (Note in the latter case, 
however, that the mass threshold falls steadily beyond 
redshift of 0.5.) 

3.7. Redshift distribution 

While a full cosmological analysis will be undertaken 
in Section 5, we briefly confirm the consistency of our 
cluster redshift distribution with expectations. As in the 
cosmological analysis, we will select our samples based on 
the signal to noise ratio of the uncorrected central decre- 
ment yo ± Syo obtained for each cluster using the filter 
corresponding to #500 = 5.'9. We first consider the S82 
clusters that have yo/Syo > 4, over the redshift range 
0.2 < z < 0.8. Secondly we consider the "cosmolog- 
ical" sample of clusters, consisting of 15 clusters with 
yo/^Vo > 5.1 and z > 0.2. The cumulative number den- 
sity as a function of redshift is shown in Figure 12. For 




0.6 0.8 

Redshift 

Figure 11. Estimate of the mass (Msooc) above which the ACT 
cluster sample within S82 is 90% complete (see Section 3.6). Lower 
panel assumes a UPP-based scaling relation with 20% intrinsic 
scatter; solid line is for S/N > 4 (full S82 sample, valid to z < 0.8), 
dotted line is for the S/N > 5.1 subsample (valid to2K 1.4). The 
upper panel shows analogous limits, but assuming scaling rela- 
tion parameters obtained for the B12 model (Section 3.4). Circles 
(crosses) are based on filtering and analysis of B12 model clusters 
for the S/N > 4 (5.1) cut. The completeness threshold decreases 
steadily above z m 0.6 because clusters at this mass are easily re- 
solved and the total SZ signal, at constant mass, increases with 
redshift. 

each of the two subsamples, we bin the clusters into red- 
shift bins of width 0.1 and perform a maximum likelihood 
fit (assuming Poisson statistics in each bin) to estimate 
as- To facilitate comparison with cosmological results 
presented in Section 5, we assume a flat ACDM cosmol- 
ogy with Q,jm = 0.28 and n s = 0.96, and use the scaling 
relation parameters obtained for the B12 model. Cluster 
count predictions are obtained starting from the cluster 
mass function of Tinker et al. (2008), and include all se- 
lection effects (intrinsic scatter, noise, and yo/5yo cut). 
The fits yield <r$ — 0.827 for the cosmological sample, 
and as = 0.837 for the S82 sample. Both of these are 
consistent with the result of the full cosmological analy- 
sis. The S82 (respectively, cosmological) sample is well 
described by the model in the sense that 53% (61%) of 
samples drawn from the same model, and having the 
same total cluster count, are less likely. 

4. COMPARISON WITH OTHER CATALOGS 

In this section we compare the Equatorial cluster cat- 
alog and the SZ derived cluster properties to those ob- 
tained by other studies in microwave, X-ray, and optical 
wavelengths. While optical studies have good overlap 
with our sample in S82 to z < 0.6, previous X-ray and SZ 
survey data include only a small fraction of the clusters in 
our sample. We also examine the question of radio con- 
tamination of cluster decrements through a comparison 
of extrapolated fluxes near our cluster positions relative 
to random positions in the field. 

4.1. Comparison to Planck Early SZ Sample 
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Figure 12. Cumulative number counts for two subsamples of the 
full cluster catalog for which confirmation is complete. The upper 
lines are data and model counts for the S82 sample of clusters 
having ya/Syo > 4 and 0.2 < z < 0.8. The lower lines represent 
the cosmological sample of 15 clusters with fixed-scale yo/Syo > 
5.1 and z > 0.2. The model for the counts is obtained from a 
maximum likelihood fit, with only erg as a free parameter. The 
model includes a full treatment of selection effects for the sample 
under consideration. 

We compare our catalog and our derived cluster 
properties, to the catalog presented in the Planck 
all-sky Early Sunyaev-Zel'dovich cluster sample (ESZ; 
Planck Collaboration 2011a). The ESZ presents 189 
clusters, of which 4 lie within the ACT Equatorial foot- 
print, and of which 2 are detected by ACT. The two clus- 
ters detected by both Planck and ACT consist of two of 
the three clusters having ACT I500 exceeding the 50% 
completeness level of the ESZ. (The third, not matched 
to the ESZ, is RCS2 J2327.4-0204.) 

The two clusters not detected by ACT are low redshift 
clusters: Abell 2440 at z = 0.091 and Abell 119 at z 
= 0.044. Based on their integrated X-ray gas temper- 
ature measurements of 3.88 ± 0.14 and 5.62 ± 0.12 keV 
(White 2000), we estimate masses of « 4 x 10 14 and 
7 x 10 14 hj Q M@, respectively; these are well below our 
90% completeness level (Section 3.6) at these redshifts. 

For the two clusters detected by both Planck and ACT, 
a summary comparison of measured cluster properties 
may be found in Table 2. MACS J2135.2-0102 is de- 
tected by ACT, at low significance, inside S82. Abell 
2355 (ACT-CL J2135.2+0125) is detected by ACT at 
high significance (S/N = 9.3) just outside the S82 re- 
gion. The specifics of each case are discussed below. 

For ease of comparison, we convert the Planck mea- 
surement of the SZ signal within 5i?5oo through the ESZ- 
provided conversion factor Y500 = Y5R 500 /1.81. We also 
use the ESZ value for #500 (which is either determined 
from X-ray luminosity measurements, or from the SZ sig- 
nal alone) to obtain an approximate value for Afsoo c - The 
ESZ analysis makes use of the "Standard" version of the 
Universal Pressure Profile, which assumes a self-similar 
scaling relation (see the Appendix of Arnaud ct al. 2010). 
We thus re-analyze the ACT yo measurements using the 
profiles and scaling relation of the Standard UPP to es- 



timate #500, Y50Q, an d -M5ooc- The results of this analysis 
differ only slightly from the results obtained using the 
full UPP '(Table 8). 

For MACS J2135.2-0102 (z = 0.329), X-ray luminos- 
ity data was not available and the ESZ presents the an- 
gular scale of the cluster based on SZ data alone. The 
scale, #500 = 1-6 ± 1.0' is the smallest #500 in the ESZ, 
and corresponds to a very low mass («lx 10 14 ). Such 
a mass seems inconsistent with the SZ signal observed 
by either Planck or ACT. Lensed submillimetcr galax- 
ies have been observed near this cluster (Ivison et al. 
2010), and the ESZ notes include a reference to possi- 
ble point source contamination. The ACT measurement 
is likely to be less contaminated by such emission, since 
we use only 148 GHz data where dusty sources are com- 
paratively dim. Overall this cluster is a peculiar case, 
and it is difficult to draw any useful conclusions from the 
disagreement of ACT and Planck measurements without 
more detailed X-ray information. 

Abell 2355 (ACT-CL J2135.2+0125) is one of the 
most significant ACT detections, and one for which the 
ACT analysis implies a very high mass. A spectro- 
scopic redshift of 0.1244 for this cluster has been ob- 
tained by Kowalski et al. (1983), and has subsequently 
been adopted in both the MCXC and the Planck ESZ. 
However, Sarazin et al. (1982) identify z = 0.231 as a 
more probable spectroscopic redshift, and this value is 
adopted by Menanteau et al. (2012b), who find it to be 
much more consistent with their photometric estimate 
of z = 0.25 ± 0.01. Furthermore, the NED 4 entry for 
this cluster refers to an unpublished spectroscopic red- 
shift z = 0.228 obtained from 3 galaxies. In order to 
compare the Planck and ACT measurements, we correct 
the Planck SZ measurements to z = 0.231. 

The cluster angular scale used by Planck is obtained 
from X-ray luminosity based masses in the MCXC. We 
compute a new mass estimate using the MCXC scaling 
relations, correcting the X-ray luminosity for the changes 
in luminosity distance and K-correction (according to the 
T = 5 keV tabulation of Bohringer et al. 2004). The 
resulting inferred mass is more than double the estimate 
at z = 0.1244. The corresponding #500 is slightly smaller, 
and so we obtain a crude estimate of the Y500 that Planck 
might have measured if they had used this angular scale. 
From the inspection of Figure 9 of Planck Collaboration 
(2011a), the axis of degeneracy for the scale and signal 
measurements lies along Y500 oc 6^ 00 with a in the range 
of 0.75 (for resolved clusters) to 1.5 (unresolved clusters). 
We compute the new value assuming a = 1, and add 
20% error to account for the uncertainty in a. This gives 
y 500 = (15 ± 4) x 10~ 4 arcmin 2 . The ACT mass and 
Y500 are in good agreement with the X-ray mass and our 
estimate of the resulting Planck SZ signal. 

Comparison to SZ measurements of three more clusters 
detected by Planck may be found in the next section. 

4.2. Comparison to Weak Lensing Masses 

In this section we examine weak lensing mass mea- 
surements of 4 clusters in the ACT Equatorial sample. 
While one is a high redshift cluster discovered by ACT, 
the other three are well-known moderate redshift clusters 

4 NASA Extragalactic Database; http://ned.ipac.caltech.edu/. 
Retrieved July 15, 2012. 
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Table 2 

Comparison of cluster properties as determined by Planck and by ACT, for two clusters from the ESZ (Planck Collaboration 2011a, 
Section 4.1) and three clusters from (Planck Collaboration 2012a, Section 4.2). For Planck, 6500 is derived from X-ray measurements of M500 
for all clusters except MACS J2135.2— 0102, for which the angular scale was determined from the SZ data only. Planck values for Abell 2355 
are corrected to redshift 0.231 as discussed in the text. For the ACT measurements, the angular scale, mass and Y500 of each cluster are 
obtained simultaneously from the 148 GHz data assuming that the cluster pressure profile is described by the Standard UPP. 



Cluster ID 


Redshift 




Planck 






ACT 








#500 


y 500 


M500 


#500 


V500 


Aisoo 






(arcmin) 


(10 -4 arcmin 2 ) 


(10 14 h£ Mq) 


(arcmin) 


(10~ 4 arcmin 2 ) 


(10 14 Mq) 


Early SZ clusters 
















MACS J2135. 2-0102 


0.325 


1.6 ± 1.0 


9.8 ± 1.9 


0.8 ± 0.8 


3.1 ± 0.4 


2.5 ± 1.2 


2.7 ± 1.0 


Abell 2355 


0.231 


5.1 


15 ± 4 


5.2 


5.3 ± 0.2 


14.3 ± 2.4 


6.3 ± 1.3 



Planck Collaboration (2012a) 

Abell 267 0.235 5.5 ± 0.3 9.3 ± 2.3 6.3 ± 1.0 5.4 ± 0.2 13.1 ± 2.4 5.6 ± 1.2 

RXC J2129.6+0005 0.234 4.7 ± 0.2 7.8 ± 2.0 4.3 ± 0.5 5.2 ± 0.2 11.4 ± 2.5 5.2 ± 1.2 

Abell 2631 0.275 5.4 ± 0.6 15.5 ± 2.3 9.8 ± 3.3 4.8 ± 0.2 11.5 ± 2.2 6.0 ± 1.3 



that have been observed by XMM-Newton and Planck. 

Weak lensing measurements of ACT-CL J0022.2-0036 
(z = 0.81) are presented in Miyatake ct al. (2012). Sub- 
aru imaging is analyzed and radial profiles of tangential 
shear are fit with an NFW profile. They obtain a mass 
estimate of M 50 o c = 8A±f% x 10 14 h^M e . While consis- 
tent with our SZ masses for any of the three model scaling 
relations, this mass is higher than the one deduced from 
the UPP scaling relation parameters and more consis- 
tent with the B12 and Nonthermal20 models. The SZ 
and lensing mass estimates are also consistent with the 
dynamical mass reported for ACT-CL J0022. 2-0036 in 
Menanteau et al. (2012b). 

The ACT Equatorial sample also includes 3 clus- 
ters (A267, A2361, and RXCJ2129.6+0005) treated by 
Planck Collaboration (2012a, hereafter PI3) in a compar- 
ison of weak lensing mass, X-ray mass proxies, and SZ 
signal. X-ray data are obtained from the XMM-Newton 
archive, and weak lensing masses for the clusters we con- 
sider here originate in Okabe et al. (2010). Each of these 
clusters is detected by ACT inside S82 with S/N > 8. 

Comparing the ACT UPP based SZ masses to the weak 
lensing masses for these three clusters we obtain a mean 
weighted mass ratio of M^/Mwl = 1.3 ± 0.2. This is 
consistent with the mass ratio found by PI3 between X- 
ray masses and weak lensing masses for their full sample 
of 17 objects. For the B12 scaling relation parameters 
the ratio is M|^/M WL = 1.8 ± 0.3. 

PI3 also report Y500, measured from Planck's multi- 
band data using a matched filter, with the scale (#500) of 
the cluster template fixed using either the X-ray or the 
weak lensing mass. For each object we see agreement in 
the Y500 measurements at the 1 to 2 — a level, and our 
mean ratio is consistent with unity. For #500 determined 
from X-ray (weak lensing) mass we find weighted mean 
ratio YUPP CT /Y bm , Planck = 0.90 ± 0.16 (0.98 ± 0.16). 

The angular scales, Y500 and M500C obtained by Planck 
and ACT are provided in Table 2. The ACT values are 
computed for the Standard version of the UPP, following 
our treatment of the ESZ clusters in section 4.1; these 
results are almost indistinguishable from those obtained 
with the full UPP treatment. 

4.3. Comparison to SZA Measurements 

Higher resolution SZ data can be obtained through in- 
terfcrometric observations. Reese et al. (2012) present 



Sunyaev-Zel'dovich Array (SZA) observations of two 
ACT Equatorial clusters at 30 GHz. For the high red- 
shift, newly discovered cluster ACT-CL J0022-0036, a 
GNFW profile is fit to the 30 GHz SZ signal, and is 
used to infer the cluster mass assuming an NFW density 
profile and virialization of the gas. This yields a mass 
estimate of Af 500c = 7.3±J;g x 10 14 h^M Q . For Abell 
2631 (ACT-CL J2337.6+0016), X-ray and SZ data are 
both used to constrain the density profile, producing a 
mass estimate of M 50 oc = 9-4±|;| x 10 14 h^M Q . These 
masses are somewhat higher than the ACT results for the 
UPP scaling relation, and are more consistent with the 
masses arising from the B12 scaling relation parameters. 

4.4. Optical Cluster Catalogs 

The extensive overlap of the ACT observations de- 
scribed in this work with the SDSS means that there 
are a number of existing optically selected cluster cata- 
logs with which the ACT SZ selected cluster sample can 
be compared (see also Menanteau et al. 2012b). For all 
the comparisons described below, we matched each cata- 
log to the ACT cluster sample using a 0.5 Mpc matching 
radius, evaluated at the ACT cluster redshift. 

Several optical cluster catalogues have been extracted 
from the SDSS legacy survey (e.g., Goto et al. 2002; 
Miller et al. 2005; Koester et al. 2007; Szabo et al. 2011; 
Wen et al. 2009, 2012). For the purposes of this compar- 
ison, we focus on the MaxBCG catalog (Koester et al. 
2007) and its successor the GMBCG catalog (Hao et al. 
2010). Both of these catalogs make use of the color- 
magnitude red-sequence characteristic of the cluster 
early type galaxy population, plus the presence of a 
Brightest Cluster Galaxy (BCG), to identify clusters. 
The MaxBCG catalog contains 13,823 clusters, of which 
492 fall within the footprint of the 148 GHz map used 
in this work, and is thought to be > 90% complete and 
> 90% pure for clusters with N ga \ > 20 over its entire 
redshift range (0.1 < z < 0.3). The GMBCG catalog 
builds on this work using the entire SDSS DR7 survey 
area, and is thought to have > 95% completeness and pu- 
rity for clusters with richness > 20 galaxies and z < 0.48. 
A total of 1903 of the 55,424 GMBCG clusters fall within 
the ACT footprint. 

We find that the ACT cluster catalog contains 8 clus- 
ters in common with MaxBCG, and 16 clusters in com- 
mon with the GMBCG catalog. There are no ACT clus- 
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ters at z < 0.3 in the SDSS DR7 footprint that are not 
cross matched with MaxBCG objects; however, there 
are 2 ACT clusters at z < 0.48 (ACT-CL J0348.6-0028 
and ACT-CL J0230.9-0024) that arc not cross matched 
with objects with richness > 20 in the GMBCG cat- 
alog within the common area 5 between the two sur- 
veys. ACT-CL J0348.6-0028 (z = 0.29) is an opti- 
cally rich system (N ga i = 56.9 ± 7.5; as measured by 
Menanteauet al. 2012b), while ACT-CL J0230.9-0024 
is optically fairly poor (N ga j = 19.9 ± 4.5). In both 
cases, Menanteau et al. (2012b) find the BCG to have 
a small offset from the SZ position (0.1 Mpc for J0348.6 
and 0.16 Mpc for J0230.9). However, neither of these 
objects has a plausible cross-match in the full GMBCG 
catalog, although we note that there is a GMBCG cluster 
(J057.14850-00.43348) at z = 0.31 located within a pro- 
jected distance of 0.65 Mpc of ACT-CL J0348.6-0028. 

We also compared the ACT catalog with that of 
Geach, Murphy, & Bower (2011). This catalog is con- 
structed using the Overdense Red-sequence Cluster Al- 
gorithm (ORCA; Murphy, Geach, & Bower 2012) and is 
the first optical cluster catalog available based on the 
deep (r re 23.5 mag) SDSS S82 region. It reaches to 
higher redshift (z re 0.6) than the catalogs based on 
the SDSS legacy survey data (such as GMBCG), and 
all of the MaxBCG clusters within a 7dcg 2 test area 
are re-detected (Murphy et al. 2012). We find that 26 
z < 0.6 ACT clusters are cross-matched with objects 
in the Geach et al. catalog. However, there are 24 ACT 
clusters, which we have optically confirmed using the S82 
data, that were not detected by Geach et al., and most 
of these (19 objects) are at z > 0.6. We also confirmed a 
further 4 objects at z > 1 in the S82 region with the ad- 
dition of K s -band imaging obtained at the Apache Point 
Observatory. This suggests that the ACT cluster catalog 
has a higher level of completeness for massive clusters at 
high redshift, although we note that the efficiency of the 
optical cluster detection algorithm used by Geach et al. 
at such redshifts is limited by the lack of infrared data. 

4.5. X-ray Cluster Catalogs 

Clusters are detected in X-rays through thermal 
bremsstrahlung emission from their hot gas atmospheres, 
and so X-ray selected cluster surveys are complemen- 
tary to SZ searches. However, current large area X- 
ray surveys are relatively shallow. The REFLEX clus- 
ter catalog (Bohringer et al. 2004) is derived from the 
ROSAT All Sky Survey data (Voges et al. 1999) and 
overlaps completely with the ACT survey. The full cata- 
log contains a total of 448 clusters reaching to z = 0.45, 
from which 17 objects fall within the footprint of the 
ACT Equatorial maps. We find that only five of these 
clusters are cross-matched with ACT cluster detections. 
The detected objects are all luminous systems, with 
Lx > 3.3 x 10 44 ergs _1 . The undetected objects are all 
lower luminosity (hence lower mass) and are at z < 0.1, 
where SZ completeness is low. 

As discussed in Section 4.1, the Planck ESZ relies on 
the M500C values presented in the MCXC (derived from 

5 To determine the overlap between the ACT maps and various 
SDSS data releases, we make use of the angular selection function 
from http: //space .mit . edu/~molly/mangle/download/ data. html 
(see Swanson et al. 2008; Hamilton & Tcgmark 2004; 
Blanton ct al. 2005, for details). 
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Figure 13. Ratio of ACT SZ determined masses to X-ray lu- 
minosity based masses from the MCXC. ACT masses assume the 
UPP scaling relation parameters. Error bars on mass include un- 
certainty from the ACT SZ measurements, and 50% uncertainty 
on MCXC masses. The weighted mean ratio is 1.03 ± 0.19 for the 
Equatorial clusters and 0.83 ± 0.13 for the full sample. 

the Lx — M relation in Arnaud et al. (2010), which has 
intrinsic scatter re 50%) to constrain the angular scale of 
detected clusters, and reduce the uncertainty in Ysr 500 . 
The MCXC catalog includes 9 clusters from the ACT 
Equatorial cluster sample, and 6 clusters from the ACT 
Southern cluster sample. 

Comparing the masses we derive from the UPP scal- 
ing relation to the MCXC masses, we obtain a mean ratio 



M so6c ■ / M 500?° = 1.03±0.19. Adding in the South- 
ern clusters (see Appendix) we find a ratio of 0.83 ±0.13. 

While the Equatorial cluster result is consistent with 
the UPP scaling relation, the full sample prefers higher 
masses. For the B12 scaling relation parameters and the 
full sample we find M^ J' B12 /M^ C = 1 u ± .17. 
For the Equatorial sample alone, this ratio is 1.41 ±0.27. 
The mass ratios are shown in Figure 13. 



4.6. Radio Point Sources 

Relative to the field, galaxy clusters are observed 
to contain an excess of radio point sources (e.g., 
Cooray et al. 1998). Studies of potential SZ signal con- 
tamination by radio sources have been carried out at fre- 
quencies below 50 GHz (e.g., Lin et al. 2009). There 
are, however, few such studies at 150 GHz. Recently 
Sayers et al. (2012) used Bolocam at 140 GHz to study 
possible radio source contamination in 45 massive clus- 
ters. They found that the SZ signals from only 25% of 
the sample were contaminated at a level greater than 
1%. The largest contamination observed was 20% of the 
SZ signal. We study this phenomenon in the ACT clus- 
ter sample, using the FIRST catalog of flux densities 
at 1.4 GHz (White et al. 1997). The lower flux density 
limit of the catalog is at most 1 mJy, and overlaps with 
the ACT Equatorial field between R.A. of 21 h 20 m and 
3 h 20 m . 

For each source in the FIRST catalog, we extrapolate 
the flux density to 148 GHz using flux-dependent spec- 
tral indices computed based on a stacking analysis in 
the ACT 148 GHz and 218 GHz maps (Gralla et al., in 



ACT: SZ Selected Galaxy Clusters 



17 



prep.). For the sources of interest here, the spectral in- 
dex ranges from approximately -0.5 (at Si. 4 = lOmJy) 
to -0.8 (at Si. 4 = 100 mJy). For each of the 63 galaxy 
clusters lying in the FIRST survey area, we take all 
sources within 2' of the cluster position and sum the 
predicted flux at 148 GHz. We find a mean flux density 
of S148 = 0.94 mJy, and the 9 most potentially contami- 
nated clusters have flux density between 2.5 and 4.2 mJy. 
Converting these to a peak brightness in CMB temper- 
ature units using the 148 GHz beam, we obtain contam- 
ination of 11 u.K on average and 50 u.K at worst. These 
levels are substantially higher than would be expected 
if total radio flux were not correlated with cluster posi- 
tion. For random positions, the contamination is only 
S148 = 0.2 mJy on average, and 1% of locations have 
total flux exceeding 2.4 mJy. 

For the purposes of inferring masses, we are interested 
in the impact the point sources may have in the measure- 
ment of the uncorrected central decrement. To explore 
this, we create a simulated ACT map of the sources us- 
ing their inferred flux densities at 148 GHz and the ACT 
beam shape. This map is then filtered with the same 
matched filter ^ described in Section 3.1. We obtain a 
prediction for the contamination level of each cluster by 
taking the maximum value in the filtered map within 2' of 
the cluster position. More than half of the clusters have 
predicted contamination less than 13 u.K and 90% have 
less than 30 uK (which is smaller than the typical mea- 
surement uncertainty). The worst contamination predic- 
tion is associated with ACT-CL J0104. 8+0002, at a level 
of 59 u.K (corresponding to 35% of the signal strength, 
and larger than the cluster noise level of 41 |xK). 

Since such contamination is difficult to model in de- 
tail without knowing the spectral indices of individual 
sources, and since the rate of significant contamination 
seems to be quite low, we make no attempt to correct for 
this effect. This may introduce a small, redshift depen- 
dent bias into scaling relation parameters obtained from 
these SZ data. 

5. COSMOLOGICAL INTERPRETATION 

Cluster count statistics, such as the number density of 
clusters above some limiting mass, are particularly sensi- 
tive to the total matter density (f2 m ) and the amplitude 
of density fluctuations (as parametrized by, e.g., erg). Be- 
cause SZ selected cluster samples can reach to arbitrar- 
ily high redshift, they also probe parameters, such as the 
dark energy equation of state parameter w, that describe 
the recent expansion history. 

In order to constrain cosmological parameters, we in- 
corporate our sample into a Bayesian analysis and com- 
pute the posterior likelihood of cosmological parameters 
given the cluster data. We begin by outlining the for- 
malism used to model the probability of our data given 
values of cosmological and scaling relation parameters. 
We then demonstrate the constraints achieved by com- 
bining the ACT cluster data with other data sets. 

5.1. Likelihood Formalism 

In this section we outline a formalism for determining 
the Bayesian likelihood of cosmological and scaling re- 
lation parameters given the ACT cluster measurements, 
including redshift, SZ, and dynamical mass information, 
when available. Our approach follows previous work in 



developing an expression for the probability of the clus- 
ter measurements based on the application of Poisson 
statistics to finely spaced bins in the multi-dimensional 
space of cluster observables. Such approaches naturally 
support non-trivial sample selection functions, the self- 
consistent calibration of scaling relations between ob- 
ject properties (i.e., between mass, SZ signal, and red- 
shift), and missing data (i.e., the absence of independent 
mass measurements for some detected clusters). This 
approach to the comparison of a sample of detected ob- 
jects to a number density predicted by a model is de- 
scribed by Cash (1979); Mantz ct al. (2010b) present a 
useful general formalism for dealing with cluster data, 
applying it to X-ray mass and luminosity measurements 
to obtain cosmological constraints while calibrating the 
mass-luminosity scaling relation. The SZ studies of 
Sehgal et al. (2011) and SPT (e.g., Benson et al. 2011) 
have applied similar techniques to SZ and X-ray data. 
Here, we develop an approach that pays particular at- 
tention to the uncertainties in the observed quantities, 
and in which the SZ signal is interpreted through the 
PBAA approach (Section 3.1). 

We assume that the cluster data consist of a sam- 
ple of confirmed clusters, and that for each cluster an 
uncorrected central Compton parameter yo ± Syo and 
redshift z ± 6z have been measured. The yo are ob- 
tained from maps filtered with ^g'g. In some cases, 
clusters may also have dynamical mass measurements, 
M^qq" (see Section 3.5). To parametrize these measure- 
ments in a way that is independent of cosmology, we 
define the observed dynamical mass parameter, rh = 

To compare the observed sample to predicted cluster 
number densities, we consider all clusters to possess an 
intrinsic uncorrected central Compton parameter yo' r , 
redshift z tr , and mass parameter m tr . These true intrin- 
sic quantities represent the values one would measure in 
the absence of any instrumental noise or astrophysical 
contamination (from, e.g., the CMB). For the SZ signal, 
yo tr should be thought of as the measurement we would 
make if we applied our filter ^ 5 ' 9 (i.e., the fixed-scale 
filter matched to the noise spectrum of our maps) to a 
map from which all noise and astrophysical contamina- 
tion had been removed. The true cluster mass parame- 
ter, m tr = M 50 oc h^E{z tv ) with E(z tr ) computed for the 
true cosmology, is representative of the halo mass rather 
than the mass inferred from an observational proxy (such 
as galaxy velocity dispersion). 

We proceed by obtaining the number density of galaxy 
clusters in the space of true cluster properties yo tr , m tr , 
and z tr . As described in Section 3.2, we make use of the 
cluster mass function of Tinker et al. (2008) to predict, 
for cosmological parameters 9, the number of clusters per 
unit redshift and unit mass within the area of the survey: 

n{m tr ,z tr \9) = d 2 N{< m tr , z tr ) / dz tl dm tT . (20) 

Here we use the notation n(a\/3), more commonly used 
with probability densities, to indicate the conditional dis- 
tribution of clusters with respect to variables a when 
variables f3 are held fixed. 

Given m tr , z tr , and scaling relation parameters xp = 
{A m , B, C, <7int), the distribution of Compton parameter 
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values, P(j/o r \fh tT , z tT , ip), is specified by equations (13) 
and (16). This may be used to compute the number den- 
sity of clusters in the full space of true cluster properties: 

n(y tT , m tr , z tr |0, t/>) = P(y tr |m tr , z tr , t/»)x 

n(m tr ,z tr |6»). (21) 

In order to compare our observed sample to the model, 
it is necessary to properly account for the sample selec- 
tion function, and for the effects of measurement uncer- 
tainty. This is especially important because our selection 
function depends explicitly on 5yo, and the mass func- 
tion (and thus the cluster density as a function of yo) is 
very steep. This full treatment of uncertainty allows us 
to include, in the same analysis, regions of the map that 
have quite different noise levels. 

In general, we imagine that the cluster observables 
x = (yo,fh, z,dyo,dfh,dz) are related to the true clus- 
ter properties yo r , m tr , z tr by some probability distribu- 
tion F(x|yo tr , m tr , z tr , /3 dyn ). In addition to describing 
the scatter of each variable about its true value, and 
accounting for the dynamical mass bias through the pa- 
rameter f3 dyn (which is dehned in Section 3.5), this dis- 
tribution also includes a description of what measure- 
ment uncertainties we are likely to encounter for given 
values of the true cluster properties. For example, while 
most of our sample have spectroscopic rcdshifts for which 
the measurement uncertainty is negligible, some clusters 
(particularly at high redshift) have photometric redshift 
estimates with relatively large uncertainties (6z ps 0.06). 
Although it may seem awkward to worry about mea- 
surement uncertainty for clusters that have not been de- 
tected, it is necessary, formally, to account for the dis- 
tribution of errors if the sample is defined based on ob- 
served cluster properties. In our particular case, the full 
probability distribution factors to: 

P(x|yV r ,™ tr , A/? dyn ) = P(yo|yo tr ,5y )P(5yo)x 

P(m|m tr ,5m,/3 dyn )P(<5m)x 
P{z\z tr ,Sz)P(Sz\z tr ). (22) 

The expression above encodes the following properties 
of the ACT observations: 

• For a given cluster, the measurements of yo, fh, and 
z are independent; we do not expect any covariance 
in the errors. In practice we assume cluster ob- 
servables are normally distributed about their true 
values (with the exception of fh; see next point). 

• The probability distribution of fh includes the ef- 
fects of the dynamical mass bias parameter, /3 dyn ; 
specifically the measured dynamical mass parame- 
ter fh is expected to be normally distributed about 
mean /? dyn m tr with standard deviation Sfh. 

• The distribution of yo errors, P(Sy~o), is indepen- 
dent of all true cluster properties. In practice 
P(Syo) is obtained from the histogram of the noise 
map. 

• The distribution of fh errors, P{5fh), is indepen- 
dent of true cluster properties. In practice the 
uncertainty in the dynamical masses is related to 



the number of galaxies used for the velocity disper- 
sion measurements. In any case, we do not need 
to understand this distribution in detail, because 
observed mass is not a factor in sample selection. 

• The distribution of z errors, P(Sz\z tr ), may de- 
pend on the true cluster redshift. While spec- 
troscopically measured redshifts are available for 
many sample clusters, high redshift clusters are 
more likely to have only photometric redshift es- 
timates. In practice, this distribution only enters 
when computing the prediction for the total num- 
ber of clusters observed within some volume. For 
suitably chosen sample redshift limits, the details 
will not matter (see discussion below). 

For a sample selected based on signal to noise ratio 
threshold s and observed redshift range [za 7 zb], we de- 
fine the selection function S(yo, z, 5yo) to take value unity 
when yo/<5yo > s and z G [za, Zb], and to take value zero 
otherwise. Then the predicted number density in the 
6-dimcnsional space of observables x is 

n(x|0,^/3 dyn ) = S(y o ,z,<5y o )x 

J dyo tr dm tr dz tr P(x|y tr , m tr , z tr , /3 dyn ) x 

n{$F,m*,z*\B,1>). (23) 

This cluster density function may be used to evaluate 
the extent to which the observed cluster data are consis- 
tent with the model 9,xp,f3 dyn . This is achieved, as in 
Cash (1979) by imagining a very fine binning in the space 
of observables. We take bins indexed by a centered at x a 
and having (6-dimensional) volume V a . In the limit of 
very fine bins, the total predicted counts in bin a is well- 
approximated by N a (8,tp,(3 d y D ) ps V a n(x a \6, ip, /? dyn ). 
Furthermore, the number of observed clusters in bin a, 
denoted c a , is either or 1. 

Letting D denote the set of bins in which a cluster 
has been observed (i.e., the a where c a = 1), we assume 
Poisson statistics in each bin, and obtain the probability 
of the data given the model parameters: 

P({xJ|0,^,/3 dyn )=P({c Q }|0,^,/3 dyn ) 

a 

= e -N M ,p-n J|y an(Xa) . 

(24) 

We have defined 

AT tot (0, /3 dyn ) = J2 N °< 4>> ^ yn ) 

a 

= J <fx n(x|0,V,/3 dyn ), (25) 

the total number of clusters that the model predicts will 
be detected. 

In equation (24), the product over occupied bins is only 
sensitive to the values of the density function at the lo- 
cations of the detected clusters. The volume elements 
V a depend on the data and the binning, but not on the 
cosmological or scaling relation parameters. They will 
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thus cancel exactly in any ratio of probabilities compar- 
ing different models. So we may write the likelihood of 
parameters 6,xj),{3 dyn given the cluster data {xj} as 

L{6, xj}, /3 dyn |{x 4 }) = P({xi}|0, xj,, /? dyn ) 

n"(x^,V,/? dyn ), (26) 

i 

where i indexes the clusters in the sample. 

When evaluating this expression in practice, we face 
the two related problems of computing the total clus- 
ter count prediction N tot , and of computing the number 
density n(xj|0, xj}, /3 dyn ) for each cluster in the sample. 
Certain simplifications make possible the efficient com- 
putation of these quantities. 

When computing n(xj[#, xj}, f3 dyn ), we approximate 
P(5z\z tr ) as being constant over the range of z tr under 
consideration. This is acceptable because the integral 
over z tr is restricted to the vicinity of the observed clus- 
ter redshift Zi by the distribution P(zi\z tT ,Sz). We thus 
replace P(5zi\z tr ) with p 4 = P(Szi\z tr )\ z tr =Zt . The pi 
can then be factored out of the integral over true cluster 
properties in equation (23). When evaluating the like- 
lihood in equation (26), these pi contribute a constant 
multiplicative factor that is independent of parameters 
0,xj},j3 dyn . So their contribution is irrelevant, and like 
the V a the pi may be dropped from the likelihood ex- 
pression. 

When computing A/tot i the procedure is simplified 
by first integrating over m tr to obtain the distribution 
n(yo tr , z tr \6, xj}). The integrals over m and Sfh are triv- 
ial to perform (independent of the form of P(Sm) and 
the value of /3 dyn ), because the selection function does 
not depend on fh or Sfh. We may then write N to t(d, xj}), 
dropping the dependence on f3 dyn . The integral over Syo 
can be accomplished with P(Sy~o) based on the noise map. 
Note that this is essential to properly predict the total 
cluster count in cases where the map area under consid- 
eration includes a variety of local noise levels. 

It is necessary to consider the impact of P(Sz\z tr ) in 
the evaluation of iVtot- The number of clusters within 
an observed redshift range z £ [za,zb] will be approxi- 
mately equal to the number of clusters with true redshift 
z tr G [za, zb\- The difference between these two numbers 
may be interpreted as the result of clusters "scattering" 
over the redshift boundary due to measurement uncer- 
tainty. The magnitude of this effect is related to the noise 
level Sz and the steepness of the distribution n(z) at the 
boundaries. 

Two properties of the ACT cluster sample allow us 
to avoid dealing with the details of P(Sz\z tr ). The first 
is that we have spectroscopic redshift measurements for 
all but a small number of high redshift clusters. This 
means that P(Sz\z tT ) strongly favors the case of Sz w 
at the low redshift sample boundary. Second, our upper 
redshift limit, which arises based on the depth of opti- 
cal confirmation observations, is at z = 1.4. Clusters at 
z = 1.4 are sufficiently rare that uncertainty in Sz does 
not greatly affect the total number of clusters within the 
full survey volume. For example, at z — 1.4 the pre- 
dicted cluster number density has fallen off significantly 
and the effects of redshift uncertainty can only affect the 



total predicted counts by less than 1%. This is well be- 
low the error due to sample variance. So provided we use 
a cluster sample for which SZ candidate follow-up (opti- 
cal/IR confirmation) is complete out to at least z = 1.4, 
we may disregard redshift uncertainty when computing 
iV to t, and integrate over z tr instead of z and Sz. 

For a more complicated data set (involving a large 
number of photometrically obtained redshifts) , the distri- 
bution P(Sz\z tY ) could be estimated based on the redshift 
error data in hand. 

Finally, in the case that mass data are not available 
for some or all of the clusters, we can simply integrate 
over our ignorance of fh. The model and scaling relation 
parameters give cluster density prediction, marginalized 
over the missing data, of 

n(y ,z,Sy ,Sz\e,xj}) = J d(Sffi) J dm n(x\6,xj}, (3 dyn ). 

(27) 

The likelihood expression for the mixed case is 
L{6,xP,f3 dyn \{x t }) cx e-JWCVOx 

n nffo^z^Syo^SzilO,*/}), (28) 

where M denotes the subset of clusters that have mass 
measurements . 

5.2. Parameter Constraints for Fixed Scaling Relations 

In this section we obtain cosmological parameter con- 
straints by combining the Equatorial cluster sample with 
various external data sets, with the SZ scaling relation 
parameters fixed to values indicated by the UPP pre- 
scription, and by the B12 and Nonthcrmal20 models. 
For each case, we fix the scaling relation parameters 
(A m , B, C, (Tint) to the values given in Table 1, and do not 
account for any uncertainty in these parameters. These 
results should thus be viewed as illustrating the poten- 
tial constraint achievable from these cluster data, should 
the uncertainty on the scaling relations be reduced. The 
constraints also demonstrate the sensitivity of cosmologi- 
cal parameter estimates to the different physical assump- 
tions entering each model. 

The posterior distributions of cosmological parame- 
ters are obtained through Markov Chain Monte Carlo 
(MCMC) sampling, driven by the CosmoMC software 
(Lewis & Bridle 2002), with matter power spectra and 
CMB angular power spectra computed using the CAMB 
software (Lewis et al. 2000). The ACT Equatorial data 
contribute to the likelihood according to equation (28). 
The cluster sample includes the 15 objects (identified in 
Table 7) that lie inside sample boundaries defined by 
redshift range 0.2 < z < 1.4 and signal to noise ratio 
cut yo/Syn > 5.1. As discussed in Section 2.4, the cut 
on yo/SyQ corresponds to the level above which the opti- 
cal confirmation campaign has successfully confirmed or 
falsified all SZ detections. 

Posterior distributions of parameters for data sets that 
do not include ACT cluster information arc obtained 
from chains released with the WMAP seven year re- 
sults (WMAP7; Komatsu et al. 2011). In some cases, 
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Figure 14. Constraints on ACDM cosmological parameters from 
WMAP7 (black line contours) and ACTcl+BBN+HO (without 
CMB information). Contours indicate 68 and 95% confidence re- 
gions. ACT results are shown for three scaling relations: UPP (or- 
ange contours), B12 (green lines), Nonthermal20 (violet contours). 
While any one scaling relation provides an interesting complement 
to CMB information, the results from the three different scaling 
relations span the range of parameter values allowed by WMAP 
measurements. 

chains also incorporate data from Baryon Acoustic Os- 
cillation experiments (BAO; Percival et al. 2010), and 
Type la supernovae (SNe; Hicken et al. 2009). When 
not constrained by CMB information, ACT cluster data 
are combined with results from big bang nucleosynthesis 
(BBN; Hamann et al. 2008) to constrain the baryon frac- 
tion and with Hubble constant measurements (denoted 
HO; Riess et al. 2009). The HO measurements are also 
used to restrict the parameter space of wCDM model 
studies, emphasizing the role clusters can play in such 
cases. 

We first consider a ACDM model with 7 free parame- 
ters representing the baryon and cold dark matter densi- 
ties (fib and Q c ), the angular scale of the sound horizon 
{9 a), the normalization (A s ) and spectral index (n s ) of 
the matter power spectrum, the optical depth of reion- 
ization (r), and the SZ spectral amplitude (Agz)- 

The potential for our cluster data to constrain f2 m and 
erg is demonstrated in Figure 14. In order to emphasize 
the impact of cluster studies, we have computed the pa- 
rameter likelihood for the ACT cluster data in combina- 
tion with BBN and HO only (these constrain the baryon 
density to n h h 2 = 0.022 ± 0.002 and the Hubble con- 
stant to H = 73.9 ± 3.6 km s Mpc" 1 ). We plot the 
marginalized two-dimensional distribution of parameters 
as and fi m , with contours showing the 68% and 95% con- 
fidence regions. The ACT cluster constraints, without 
CMB information, are seen to nicely complement the re- 
sults from WMAP7. But the variation in the constraints 
between models shows the degree to which uncertainty 
in the cluster physics diminishes the constraining power. 
The marginalized parameter values are provided in Ta- 
ble 3. 

In a wCDM model (which differs from the 
ACDM model in that the equation of state param- 
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Figure 15. Constraints on aiCDM cosmological parameters from 
WMAP7+H0 (solid black lines), and WMAP7+ACTcl+H0 for 
three scaling relations (B12 scaling relation is green lines; UPP 
is orange contours; Nonthermal20 is violet contours.) 

eter w may deviate from —1), cluster counts are 
sensitive to the effect of dark energy on the expansion 
rate of the recent universe. As shown in Table 4, the 
ACT cluster data provides slightly improved constraints 
on (is relative to WMAP7 alone; but the true power of 
the cluster data is to break degeneracies between as, 
O m , and w. We present composite parameters defined 
by (78(fim/0.27) - 4 and w(f2 m /0.27) to express these 
improvements . 

In Figure 15, we show 2-d marginalized constraints for 
O m ,(T8, and w. Note that, for this plot only, we have 
included the HO prior, to partially break the degeneracy 
between these three parameters and to emphasize the 
role that normalization of the SZ scaling relation plays 
in this space. 

Overall, we find no disagreement between WMAP 7 and 
the ACT cluster data for any of the three model-based 
scaling relations considered. While the three scaling re- 
lations produce results that almost completely span the 
range of D, m and as preferred by WMAP7, a better un- 
derstanding of scaling relation parameters can provide 
significant improvements in parameter constraints given 
even a relatively small cluster sample. This is addressed 
in the next section. 

5.3. Parameter Constraints from Dynamical Mass Data 

As an alternative to fixing the SZ scaling relation pa- 
rameters based on models, in this section we perform 
a cosmological analysis using the ACT Southern and 
Equatorial cluster samples, including the dynamical mass 
information for the Southern clusters from Sifon et al. 
(2012). The two samples are included as separate contri- 
butions to the likelihood. In this analysis, the dynamical 
masses directly inform the cosmology, while the selection 
function is understood in terms of the observed SZ signal, 
through the modeling of the cluster signal with the UPP. 
In this framework the scaling relation parameters will be 
naturally constrained to be consistent with the observed 
sample sizes and with the yo measurements of the clus- 
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Table 3 

Cosmological parameter constraints for the ACDM model. Numbers indicate the mean and standard deviation of the 
marginalized posterior distribution. ACTcl results for scaling relations based on the Universal Pressure Profile (UPP), and the 
B12 and Nonthermal20 (Non) models do not include marginalization over scaling relation uncertainties (Section 5.2). 
ACTcl(dyn) results use Southern and Equatorial cluster data, including dynamical mass measurements for the Southern 

clusters (Section 5.3). 
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Table 4 

Cosmological parameter constraints for the flat uiCDM model. Cluster data provide important 
constraints in the space of erg, f2 m and w. ACTcl results are presented for scaling relations based on 
the Universal Pressure Profile (UPP), and the B12 and Nonthermal20 (Non) models; these do not 
include marginalization over scaling relation uncertainties. Results constrained using dynamical mass 
data (dyn) include a full marginalization over SZ scaling relation parameters. 







Parameter (toCDM) 






Data set 


Q c h 2 


CT8 


a 8 (n m /0.27) 0A 


«j(Q m 


/0.27) 


Without ACT Cluster Data 












WMAP7 


0.111 ± 0.006 


0.832 ± 0.134 


0.790 ± 0.065 


-0.95 


± 0.13 


WMAP7 + SNe 


0.111 ± 0.006 


0.791 ± 0.042 


0.798 ± 0.060 


-0.99 


± 0.11 



Fixed Scaling Relations (§5.2) 

WMAP7 + ACTcl(UPP) 0.108 ± 0.002 0.854 ± 0.106 0.766 ± 0.018 -0.90 ± 0.06 

WMAP7 + ACTcl(B12) 0.115 ± 0.003 0.915 ± 0.111 0.849 ± 0.021 -1.04 ± 0.07 

WMAP7 + ACTcl(Non) 0.119 ± 0.003 0.915 ± 0.116 0.887 ± 0.023 -1.11 ± 0.08 

Dynamical Mass Constraints (§5.3) 

WMAP7 + ACTcl(dyn) 0.116 ± 0.005 0.921 ± 0.108 0.851 ± 0.052 -1.05 ± 0.11 

WMAP7 + ACTcl(dyn) + SNe 0.115 ± 0.004 0.835 ± 0.034 0.858 ± 0.049 -1.08 ± 0.10 



Table 5 

Cosmological parameter constraints for the flat uiCDM model, for various combinations of WMAP7, ACT cluster data 
(with scaling relation constrained using dynamical mass data), and Type la Supernovac results. 



Data set 


n c h 2 


Parameter (u>CDM) 
f2 m erg 


h 


w 


Without ACT Cluster Data 

WMAP7 

WMAP7 + SNe 
Dynamical Mass Constraints 

WMAP7 + ACTcl(dyn) 

WMAP7 + ACTcl(dyn) + SNe 


0.111 ± 0.006 
0.111 ± 0.006 

0.116 ± 0.005 
0.115 ± 0.004 


0.259 ± 0.096 
0.276 ± 0.020 

0.237 ± 0.080 
0.289 ± 0.017 


0.832 ± 0.134 
0.791 ± 0.042 

0.921 ± 0.108 
0.835 ± 0.034 


0.753 ± 0.131 
0.697 ± 0.016 

0.792 ± 0.119 
0.691 ± 0.014 


-1.117 ± 0.394 
-0.969 ± 0.054 

-1.306 ± 0.356 
-1.011 ± 0.052 



ters that also have dynamical mass measurements. 

For the Equatorial clusters we apply the same sample 
selection criteria used in Section 5.2, and thus include the 
same 15 clusters. As before, these clusters contribute to 
the likelihood through their observed redshifts and yo 
measurements . 

For the Southern cluster sample, we obtain yo measure- 
ments from the three-season 148 GHz maps as described 
in the Appendix. For the cosmological analysis we re- 
strict the Southern sample based on a signal to noise 
ratio threshold of yo/Syo > 5.7 and an observed redshift 
requirement of 0.315 < z < 1.4. The yo/Syo threshold 



is high enough to exclude new, unconfirmed candidates 
in our analysis of the three-season Southern maps. The 
lower redshift bound restricts the sample to the clusters 
for which Sifon et al. (2012) have measured dynamical 
masses. This yields a sample of seven clusters, which 
are identified in Table 9. Of the nine clusters used in 
Sehgal et al. (2011), our sample includes the five clus- 
ters at z > 0.315. The sample includes the exceptional 
cluster ACT-CL J0102— 4915; the inclusion or exclusion 
of this cluster does not change the cosmological param- 
eter constraints significantly. The contribution from the 
Southern clusters to the likelihood is in the form of equa- 
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Figure 16. Constraints on ACDM cosmological pa- 

rameters from Equatorial and Southern clusters. Re- 
sults from ACTcl(dyn)+BBN+H0 (violet contours), and 
WMAP7+ACTcl(dyn) (green contours), which both include full 
marginalization over scaling relation and dynamical mass bias 
parameters, may be compared to WMAP alone (solid black lines). 
Dotted line shows constraints for ACTcl+BBN+H0, using the 
same cluster sample but with the scaling relation fixed to the 
central values obtained from the dynamical mass fit of Section 3.5; 
note the similarity to contours in Figure 14 obtained for Equatorial 
SZ data with B12 fixed scaling relation parameters. Dashed blue 
line shows WMAP7+ACTcl(dyn), with full marginalization over 
scaling relation parameters but with fj d y n fixed to 1.33. 

tion (26). 

The data set consisting of this combination of Equato- 
rial SZ data and Southern SZ and dynamical mass data 
is denoted ACTcl(dyn). These data are combined with 
other data sets in an MCMC approach to parameter es- 
timation as described in Section 5.2, except that now we 
allow the SZ scaling relation parameters to vary, assum- 
ing fiat priors over ranges 1.7 < A m < 0.9, — 1 < B < 3, 
-2 < C < 2 and < of nt < 2. These priors are not 
intended to be informative, but rather to permit the 
scaling relation parameters to range freely over values 
supported by the data. For the dynamical mass bias 
parameter, we apply a Gaussian prior corresponding to 
l//3 dyn = 1.00 ± 0.15, 6 corresponding to the discussion 
of Evrard ct al. (2008) as described in Section 3.5. 

5.3.1. ACDM Constraints 

The effect of the increased freedom in the scaling rela- 
tion parameters may be seen in Figure 16, which shows 
the confidence regions on ag and fi m , in a ACDM model, 
from ACT cluster data combined with BBN and HO. The 
scaling relation parameters are not well constrained with- 
out some prior information, so for this chain only we 
include a Gaussian prior on the rcdshift evolution corre- 
sponding to C = 0.0 ± 0.5, based on the B12 model fits. 
Compared to results for fixed scaling relation parame- 
ters, the distribution of acceptable eg and f2 m is broader 
and skewed, at fixed fi m , towards high <jg. 

6 This description of the prior is an artifact of our initial im- 
plementation of the likelihood, where the parameter describing the 
bias corresponded to l//9 dyn . 
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Figure 17. Constraints on ACDM cosmological parameters 
from the combined Southern and Equatorial cluster samples, in- 
cluding dynamical mass measurements for the Southern clus- 
ters and full marginalization over scaling relation parameters. 
WMAP7 and WMAP7+ACTcl(dyn) are identified as in Fig- 
ure 16 (solid black line and green contours, respectively). 
Also shown are WMAP7+BAO+H0 (dotted black line) and 
WMAP7+ACTcl(dyn)+BAO+H0 (dashed blue lines). 

In a ACDM model, the addition of ACTcl(dyn) data 
improves the constraints on erg, and h relative to 
WMAP7 alone by factors of 0.8 to 0.9, as can be seen in 
Figure 16. WMAP7+ACTcl(dyn) prefers slightly larger 
values of er 8 than docs WMAP7+BAO+H0. For the 
composite parameter (Tg(fi m /0.27) 3 , the combination 
of WMAP7+BAO+H0+ACTcl(dyn) improves the un- 
certainty by a factor of 0.6 compared to WMAP7. Pa- 
rameter values are presented in Table 3 and confidence 
regions for <t§ and Q m are plotted in Figure 17. 

Within the WMAP7+ACTcl(dyn) chain for ACDM 
the dynamical mass bias parameter is pushed to /3 dyn = 
1.12 ± 0.17 (corresponding to l//3 dyn = 0.91 ± 0.12), a 
substantial change given the prior on /3 dyn . To explore 
the consequences of a large systematic bias in the dy- 
namical mass measurements, we study a ACDM chain 
run with fixed ^ dyn = 1.33. We find that cr 8 and Vl m 
move towards the central values preferred by WMAP 7, 
with parameter uncertainty slightly reduced. The confi- 
dence contours associated with this chain can be seen in 
Figure 16. 

5.3.2. Neutrino Mass Constraints 

As an extension to ACDM, we also run chains where 
the cosmic mass density of neutrinos is allowed to vary. 
Constraints are interpreted in terms of the sum of the 
neutrino mass species according to the relation Q, v h 2 = 
^^771^/(93 eV). Combining the ACT cluster data with 
WMAP7 and BAO+H0 leads to significant improve- 
ments in this constraint, as shown in Figure 18 and Ta- 
ble 6. For WMAP7+BAO+H0+ACTcl(dyn) we obtain 
an upper limit, at 95% confidence, of m v < 0.29 eV. 
The improvement in this constraint is driven by the pref- 
erence of the ACT cluster data for values of a% and 
f2 m that are in the upper range of those consistent with 
WMAP. Interpretations of this preference arc discussed 
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Figure 18. Constraints within an extension to ACDM that al- 
lows for non-zero neutrino density. The data sets shown are 
WMAP7+BAO+H0 (dotted black lines), WMAP7+ACTcl(dyn) 
(green contours), and WMAP7+ACTcl(dyn)+BAO+H0 (solid 
blue lines). The total number of rclativistic species is fixed to 
N cft = 3.046. 



below. 



5.3.3. wCDM Constraints 



Within a wCDM model, we consider the ACTcl(dyn) 
data in combination with WMAP7 and with SNe 
(which provides important, complementary con- 
straints on the recent cosmic expansion history). 
The WMAP7+ACTcl(dyn) are consistent with 
WMAP7+SNe, but with a preference (as was found in 
ACDM) for slightly higher values of a% and O m . The 
importance of cluster information is demonstrated by the 
improvement, over both WMAP7 and WMAP7+SNe, in 
the composite parameter (78(r2 m /0.27) ' 4 (see Table 4). 
As a result, the main impact of adding the cluster data 
to either WMAP7 or WMAP7+SNe is to reduce the 
uncertainties in as and f2 m by factors of « 0.8. The 
wCDM parameter constraints for each combination of 
ACTcl(dyn) and SNe with WMAP7 arc presented in 
Table 5, with marginalized 2-d confidence regions shown 
in Figure 19. 

5.3.4. Scaling Relation Constraints 

The marginalized constraints on the SZ scaling rela- 
tion parameters derived for the ACTcl(dyn) chains are 
presented in Table 1. The parameters indicating devia- 
tions from self-similar scaling with mass and redshift (B 
and C) are each consistent with 0, and consistent with 
the fits to all models. The intrinsic scatter, a lriti is some- 
what higher than the 20% seen in the model results, but 
it is not well-constrained by these data. Furthermore, 
since the data span a fairly restricted range of masses, 
there is significant covariance between A m ,B, and a- m t. 
The value of A is provided for M pivot = 7x 10 14 h^M Q , 
chosen to produce negligible covariance between A and 
B. While the cosmological results change only slightly 
when J0102— 4915 is excluded from the Southern clus- 
ter sample, we note that scaling relation parameters are 



2 -1.0 - 




0.25 0.30 



-1.2 -1.0 -0.8 0.7 0.8 0.9 1.0 

W 0"8 



Figure 19. Constraints on toCDM cosmological parameters 
from the combined Southern and Equatorial cluster samples, 
with scaling relation parameters constrained based on dynam- 
ical mass measurements. Data sets shown are WMAP7+SNe 
(dotted black lines), WMAP7+ACTcl(dyn) (green contours), and 
WMAP7+ACTcl(dyn)+SNe (solid blue lines). The units of H 
are km s — 1 Mpc . 
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Figure 20. Constraints on SZ scaling relation parameters 
from the combined Southern and Equatorial cluster sam- 
ples, ACTcl(dyn), constrained based on dynamical mass mea- 
surements in a cosmological MCMC. Creen contours arc 
for WMAP7+ACTcl(dyn) chain; black solid lines arc for 
WMAP7+ACTcl(dyn) but with J0102-4915 excluded from 
the Southern sample. Dotted line shows constraints for 
ACTcl(dyn)+BBN+H0 (i.e., without CMB information), but with 
a Gaussian prior on C of 0.0 ± 0.5. 

somewhat more affected, and in particular that the slope 
parameter B drops from 0.36 ± 0.36 to 0.06 ± 0.27 and 
the scatter er int drops from 0.42±0.19 to 0.33±0.17. The 
2-d confidence regions for scaling relation parameters are 
shown in Figure 20. 

5.3.5. Discussion 
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Table 6 

Cosmological parameter constraints for ACDM, extended with one additional parameter for non-zero 
neutrino density. The cluster data greatly assist in breaking the degeneracy between erg, f! m , and the 
neutrino density (as parametrized by J^m^). ACTcl(dyn) results use South and Equatorial cluster data, 
including dynamical mass measurements for the Southern clusters (Section 5.3). 



Parameter (ACDM + J2 m v) 
Data set Q m erg h ^2 m v ( e V) 

95% CL 



WMAP7 + BAO + HO 0.282 ± 0.018 0.742 ± 0.053 0.693 ± 0.016 < 0.58 

WMAP7 + ACTcl(dyn) 0.325 ± 0.041 0.787 ± 0.041 0.663 ± 0.029 < 0.57 

WMAP7 + ACTcl(dyn) + BAO + HO 0.289 ± 0.018 0.802 ± 0.031 0.690 ± 0.015 < 0.29 



We note that the constraints on cosmological param- 
eters are due mostly to the inclusion of the dynami- 
cal mass measurements for the Southern clusters, rather 
than due to the SZ data of the larger Equatorial sam- 
ple. However, while the removal of the Equatorial sam- 
ple from the likelihood computation produces no change 
in the cosmological parameter constraints, it leads to a 
significant weakening in the SZ scaling relation param- 
eter constraints. This is because the Equatorial sam- 
ple constrains the scaling relation parameters to values 
that predict sample selection functions consistent with 
the Equatorial sample size. By including the Equatorial 
SZ measurements in the likelihood, we obtain simultane- 
ous, self-consistent constraints on cosmological parame- 
ters and the SZ scaling relation parameters. Despite the 
low weight of the SZ information in the cosmological pa- 
rameter constraints, it remains true that the sample is 
SZ selected and thus approximately mass-limited over a 
broad range of redshifts. 

Immediate improvement in the cosmological constrain- 
ing power can be obtained through improved calibra- 
tion of the SZ mass relation. This is underway in the 
form of a campaign to collect dynamical mass data for 
the Equatorial sample (Sifon et al, in prep.). The cur- 
rent analysis would also benefit from an improved un- 
derstanding of any systematic biases in the measure- 
ment of halo mass using galaxy velocity dispersions. 
In addition to the dynamical mass data, we are pur- 
suing weak lensing (e.g., Miyatake et al. 2012), X-ray 
(e.g., Menantcau ct al. 2012a), and additional SZ (e.g., 
Reese et al. 2012) measurements to improve constraints 
on the SZ-mass scaling relation from the ACT cluster 
sample. All such mass measurements can be easily in- 
cluded in our formulation. The sample of clusters ap- 
propriate for our approach to the cosmological analysis 
will grow as targeted follow-up on the Equatorial field 
candidates is completed to lower S/N ratios. 

While this work presents a new way to quantify the 
SZ-mass relation for 10 15 M Q clusters, we note that it is 
but one component of a growing web of observation that 
tie optical and SZ data together. Recently the Planck 
team presented results (Planck Collaboration 2012c) on 
the SZ emission from SDSS galaxies that extended the 
SZ-mass relation down to 10 13 M Q systems and showed 
that the gas properties of dark matter haloes are similar 
to those in massive clusters. 

6. CONCLUSION 

The ACT Equatorial maps at 148 GHz, along with 
an optical and infrared confirmation campaign, have 
yielded a sample of 68 confirmed, SZ selected galaxy 
clusters in the rcdshift range from approximately 0.1 to 



1.4, in an area of 504 square degrees. Inside the 270 
square degree overlap with SDSS Stripe 82, the sam- 
ple is estimated to be 90% complete above a mass of 
M^qoc pa 4.5 x 10 14 H^qMq (though this number is de- 
pendent on the normalization of the SZ-mass scaling re- 
lation) . 

In order to use the SZ signal to determine cluster mass 
and to constrain cosmological parameters, we have de- 
veloped a framework that predicts the cluster amplitude 
in a matched-filtered map based on cluster mass and red- 
shift. The approach naturally handles varying noise lev- 
els in a map, and can be adapted easily to accommodate 
alternative cluster pressure profiles and scaling relations. 

The framework is based on the Universal Pressure Pro- 
file of A10, but we obtain alternative normalizations 
through an analysis of several cluster physics models. 
While the normalization obtained directly from the UPP 
is consistent with X-ray masses from the MCXC, the 
calibration to the B12 model is more consistent with dy- 
namical mass measurements of the ACT Southern clus- 
ter sample. From this we conclude that it is necessary 
to consider a broad range of scaling relation parameters 
when using the cluster data to constrain cosmological 
parameters. 

We have presented constraints in ACDM and wCDM 
models using the ACT Equatorial cluster sample for sev- 
eral fixed normalizations of the SZ scaling relations. In 
each case the cosmological constraints are consistent with 
WMAP. The results obtained for fixed scaling relations 
demonstrate the potential for such a cluster sample to 
provide important new cosmological information, even 
for modest cluster samples, provided the systematics of 
the SZ-mass calibration can be better understood. 

We have also demonstrated cosmological constraints 
based on the combination of SZ measurements and dy- 
namical mass data, in which the four parameters of the 
SZ scaling relation are calibrated simultaneously with 
cosmological parameters. The results provide significant 
constraints that are complementary to CMB, BAO, and 
Type la Supernovae data. The scaling relation parame- 
ters obtained in this analysis are consistent with the B12 
model and inconsistent with the normalization arising 
directly from the UPP. 

Our results are consistent with a study of the SZ sig- 
nal from the ACT Southern cluster sample (Schgal ct al. 
2011), as well as the ag constraints from the skew- 
ness analysis of the ACT Equatorial 148 GHz maps. 
(Wilson et al. 2012). 

While the ACT SZ cluster data prefer matter density 
parameters that are at the upper end of those supported 
by WMAP7 data, our results are consistent with cluster 
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studies that incorporate a variety of other mass proxies. 
They are also in agreement with the recent ACT angular 
power spectrum results (Sievers et al, in prep.), which 
also favor the upper limits of as and fi m permitted by 
WMAP7. 

The combined X-ray, / g as 7 Supcrnovae, BAO, and 
CMB results of Mantz et al. (2010b) produce cr 8 = 0.80± 
0.02 (0.79±0.03) and Q m = 0.257±0.015 (0.272±0.016) 
in ACDM (iuCDM); these results are marginally consis- 
tent with but lower than the results of our analysis. For 
the ACTcl(dyn)+HO+BBN run in ACDM, we compute 
cr 8 (fi m /0.25) - 47 = 0.87 ±0.04, which is higher than, but 
consistent with, the value of 0.813 ± 0.027 obtained for 
the X-ray cluster study of Vikhlinin et al. (2009b). 

Using SPT cluster SZ and X-ray mass information, 
Reichardt et al. (2012b) obtain as and fl m measurements 
in a combined WMAP7, high-£ power spectrum, and 
clusters analysis that lie slightly below the central values 
preferred by CMB measurements alone (including Siev- 
ers et al, in prep.). From their Figure 5 we also estimate 
that their cluster data in combination with HO and BBN 
(i.e., without CMB information) produce a composite pa- 
rameter constraint of <78(£l m /0.27) 03 = 0.77±0.05, lower 
by roughly 1.3-cr than either our WMAP7+ACTcldyn or 
BBN+HO+ACTcldyn results. Our results are in general 
agreement with SPT despite several differences in the 
SZ signal interpretation and mass calibration, which we 
summarize here. 

SPT makes use of cluster signal simulations, analo- 
gous to the B12 models used here, in order to interpret 
an observable based on signal to noise ratio. In con- 
trast, our analysis relies instead on the assumption of 
a simple relation between mass and the cluster pressure 
profile. We test the PBAA approach on simulated maps 
based on models of cluster physics, but we do not use 
these models to place priors on the scaling relation pa- 
rameters. The SPT mass calibration is ultimately de- 
rived from measurements of Yx, defined as the product 
of the core-excised X-ray temperature and the cluster 
gas mass (Kravtsov et al. 2006). The use of Yx as a 
mass proxy has, in turn, been calibrated to weak lensing 
masses (Vikhlinin et al. 2009a). The slight preference of 
our data for higher values of as and fl m (for both the 
ACTcl(dyn) and ACTcl(B12) studies) compared to SPT 
may indicate a complicated relationship between the var- 
ious mass proxies. 
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Table 7 

Table of confirmed galaxy clusters in the ACT Equatorial region. Coordinates 

are in the J2000 standard equinox. Rcdshifts for which uncertainties are 
quoted are photometric; all others are spectroscopic (full details can be found 
in Mcnantcau et al. 2012b). The Region (Reg.) column indicates whether the 
cluster lies within the coverage of SDSS Stripe 82 or in the shallower coverage 
of SDSS DR8; an asterisk denotes clusters used for the cosmological analysis 
of Section 5. The signal to noise ratio and filter scale 6500 are provided for 
the matched filter template that yielded the most significant detection. The 
uncorrected central Compton parameter yo is obtained from maps filtered 
with the matched filter having 8500 = 5'9 (see Section 3.1). 
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Abell 2356 (3) 

GMBH J035.40587-00.21967 (5) 
GMBH J335. 19871-00. 69024 (5) 

RCS2 J2327.4-0204 (8) 
Abell 2355 (3) 
Abell 370 (3) 

WHL J205853.1+012411 (1) 
WHL J004512. 5-015232 (1) 
RXC J2050.7+0123 (9) 
WHL J212823.4+013536 (1) 
WHL J202512.8+003134 (1) 

WHL J230739.9+013056 (1) 
Abell 2397 (3) 
RXC J0301.6+0155 (9) 
RXC J2051. 1+0216 (9) 
Abell 409 (3) 

WHL J015624.3-012317 (1) 
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ACT ID 


R.A. 


Dec. 


Redshift 


Reg. 


S/N 


#500 


yo 


Alternate ID (ref.) 






(°) 


(°) 








(arcmin) 


(io- 4 ) 




ACT- 


-CL J0219.9+0129 


34.9759 


1.4973 


0.35 ± 0.02 


DR8 


4.9 


2.35 


0.62 ± 0.14 


NSCS J021954+013102 (10) 


ACT- 


-CL J0240. 0+0116 


40.0102 


1.2693 


0.62 ± 0.03 


DR8 


1.8 


4.71 


0.70 ± 0.15 


WHL J024001. 7+011606 (1) 


ACT- 


-CL J0008. 1+0201 


2.0418 


2.0204 


0.36 ± 0.04 


DR8 


4.7 


4.71 


0.96 ± 0.21 


WHL J000810.4+020112 (1) 


ACT- 


-CL J0139.3-0128 


24.8407 


-1.4769 


0.70 ± 0.03 


DR8 


4.3 


1.18 


0.54 ± 0.17 





Note. — References: (1) Wen et al. (2012); (2) Hao et al. (2010); (3) Abell (1958); (4) Ebcling et al. (1998); (5) Geach et al. (2011); 
(6) Wen et al. (2009); (7) Goto et al. (2002); (8) Gralla ct al. (2011); (9) Bohringer et al. (2000); (10) Lopes et al. (2004). 



Table 8 

SZ-derived mass estimates for ACT clusters, inferred as described in 
Section 3.2. The yo-M scaling relation is fixed to either the UPP result, the 
fit to the B12 model, the fit to Nonthcrmal20 model (superscript "non"), or 
the fit to dynamical masses of Sifon ct al. (2012) (superscript "dyn"). Masses 
and SZ quantities are computed from the cluster's uncorrected central 
Compton parameter y~o and redshift measurements, including correction for 
mass function bias. The cluster scale 6500 , value of the bias function Q(m, z), 
and integrated Compton parameter Y500 are those inferred from the UPP 

scaling relation. Uncertainties do not include uncertainty in the scaling 
relation parameters, but do include the effects of intrinsic scatter and the 
measurement uncertainty. SZ and mass errors are highly correlated. 





ID 


#500 
(arcmin) 


Q( 


in. 


*) 


(10- 


*500 
arcmin ) 


(10 


JU 500c 


> f B12 

JU 50Qc 

(W 14 h^M Q ) 




M dyn 
(lO^h^Mo) 


ACT-CL 


J0022. 2-0036 


2.1 


± 


0.1 


0.52 


± 


0.02 


3.8 


± 0.6 


5.5 


± 





,9 


7.3 


± 1.2 


8.0 


± 


1.4 


7.7 


± 1.6 


ACT-CL 


J0326.8-0043 


3.1 


± 


0.1 


0.75 


± 


0.02 


5.5 


± 0.9 


5.5 


± 


1 





7.4 


± 1.4 


8.9 


± 


1.6 


7.9 


± 1.9 


ACT-CL 


J0152. 7+0100 


5.4 


± 


0.2 


0.98 


± 


0.01 


13.0 


± 2.3 


5.7 


± 


1 


1 


7.9 


± 1.6 


9.7 


± 


1.9 


8.7 


± 2.3 


ACT-CL 


J0059.1-0049 


2.1 


± 


0.1 


0.53 


± 


0.02 


3.5 


± 0.6 


5.2 


± 





.9 


6.9 


± 1.2 


7.6 


± 


1.3 


7.2 


± 1.5 


ACT-CL 


J2337.6+0016 


■1.8 


± 


0.2 


0.94 


± 


0.01 


11.5 


± 2.2 


6.1 


± 


1 


.2 


8.1 


± 1.7 


10.3 


± 


2.0 


9.3 


± 2.5 


ACT-CL 


J2129.6+0005 


5.2 


± 


0.2 


0.97 


± 


0.01 


11.4 


± 2.4 


5.3 


± 


1 


.1 


7.3 


± 1.6 


9.1 


± 


1.9 


7.9 


± 2.2 


ACT-CL 


J0014.9-0057 


2.8 


± 


0.1 


0.69 


± 


0.02 


5.0 


± 0.9 


5.7 


± 


1 


1 


7.6 


± 1.4 


9.0 


± 


1.7 


8.2 


± 1.9 


ACT-CL 


J0206.2-0114 


2.2 


± 


0.1 


0.54 


± 


0.02 


2.6 


± 0.6 


1.3 


± 





,8 


5.7 


± 1.1 


6,1 


± 


1.2 


5.7 


± 1.4 


ACT-CL 


J0342.0+0105 


1.5 


± 


0.1 


0.37 


± 


0.02 


1.8 


± 0.4 


3.7 


± 





7 


1.8 


± 1.0 


4.8 


± 


1.0 


1.5 


± 1.1 


ACT-CL 


J2154.5-0049 


2.7 


± 


0.2 


0.68 


± 


0.03 


3.4 


± 0.9 


4.3 


± 





.9 


5.7 


± 1.3 


6.8 


± 


1.5 


5.7 


± 1.6 


ACT-CL 


J0218.2-0041 


2.1 


± 


0.1 


0.54 


± 


0.03 


2.2 


± 0.6 


3.8 


± 





.8 


4.9 


± 1.1 


5.6 


± 


1.3 


1.7 


± 1.3 


ACT-CL 


J0223.1-0056 


2.1 


± 


0.1 


0.53 


± 


0.03 


2.2 


± 0.6 


3.8 


± 





.8 


5.0 


± 1.1 


5.7 


± 


1.3 


1.8 


± 1.4 


ACT-CL 


J2050.5-0055 


2.2 


± 


0.1 


0.56 


± 


0.04 


2.2 


± 0.7 


3.8 


± 





9 


1.9 


± 1.2 


5.6 


± 


1. 1 


1.7 


± 1.4 


ACT-CL 


J0044.4+0113 


1.3 


± 


0.1 


0.32 


± 


0.04 


1.1 


± 0.4 


2.7 


± 





.8 


3.5 


± 1.0 


3.5 


± 


1.1 


3.0 


± 1.1 


ACT-CL 


J0215.4+0030 


1.9 


± 


0.1 


0.47 


± 


0.03 


1.8 


± 0.5 


3.5 


± 





.8 


4.5 


± 1.1 


5.0 


± 


1.2 


1.2 


± 1.3 


ACT-CL 


J0256. 5+0006 


3.1 


± 


0.2 


0.76 


± 


0.03 


3.4 


± 1.0 


3.8 


± 





.9 


5.0 


± 1.2 


6.1 


± 


1.4 


1.8 


± 1.5 


ACT-CL 


J0012.0-0046 


1.2 


± 


0.1 


0.28 


± 


0.03 


1.1 


± 0.5 


3.0 


± 





.8 


3.9 


± 1.1 


3.7 


± 


1.0 


3.3 


± 1.2 


ACT-CL 


J0241. 2-0018 


1.8 


± 


0.2 


0.44 


± 


0.06 


1.1 


± 0.5 


2.5 


± 





.9 


3.2 


± 1.2 


3.6 


± 


1.3 


2.7 


± 1.2 


ACT-CL 


J0127.2+0020 


3.0 


± 


0.2 


0.73 


± 


0.04 


2.8 


± 0.9 


3.3 


± 





.9 


■1.1 


± 1.2 


5.3 


± 


1.4 


1.1 


± 1.4 


ACT-CL 


J0348.6+0029 


3.5 


± 


0.3 


0.81 


± 


0.04 


3.3 


± 1.2 


3.1 


± 





.9 


4.1 


± 1.2 


5.0 


± 


1.5 


3.7 


± 1.4 


ACT-CL 


J0119.9+0055 


1.9 


± 


0.1 


0.47 


± 


0.04 


1.7 


± 0.5 


3.3 


± 





,8 


1.3 


± 1.1 


1.7 


± 


1.3 


3.9 


± 1.3 


ACT-CL 


J0058.0+0030 


1.8 


± 


0.1 


0.45 


± 


0.04 


1.5 


± 0.5 


3.2 


± 





,8 


1.1 


± 1.1 


1.5 


± 


1.3 


3.7 


± 1.3 


ACT-CL 


J0320.4+0032 


3.0 


± 


0.2 


0.72 


± 


0.04 


2.8 


± 0.9 


3.3 


± 





.9 


1.1 


± 1.2 


5.3 


± 


1.1 


1.1 


± 1.4 


ACT-CL 


J2302. 5+0002 


2.5 


± 


0.2 


0.62 


± 


0.04 


2.5 


± 0.8 


3.7 


± 





.9 


1.8 


± 1.3 


5.6 


± 


1.5 


4.6 


± 1.5 


ACT-CL 


J2055.4+0105 


2.9 


± 


0.2 


0.70 


± 


0.04 


2.8 


± 1.0 


3.5 


± 





.9 


1.6 


± 1.3 


5.5 


± 


1.5 


1.3 


± 1.5 


ACT-CL 


J0308.1+0103 


1.9 


± 


0.2 


0.48 


± 


0.05 


1.3 


± 0.6 


2.7 


± 





.8 


3.1 


± 1.1 


3.8 


± 


1.3 


3.0 


± 1.3 


ACT-CL 


J0336.9-0110 


1.2 


± 


0.1 


0.27 


± 


0.03 


1.0 


± 0.4 


2.5 


± 





7 


3.2 


± 0.9 


3.1 


± 


0.9 


2.7 


± 1.0 


ACT-CL 


J0219.8+0022 


2.2 


± 


0.2 


0.57 


± 


0.05 


1.7 


± 0.7 


3.0 


± 





.9 


3.9 


± 1.2 


1.5 


± 


1.4 


3.5 


± 1.4 


ACT-CL 


J0348.6-0028 


3.1 


± 


0.3 


0.75 


± 


0.05 


2.6 


± 1.0 


3.0 


± 





.9 


1.0 


± 1.2 


4.8 


± 


1.5 


3.5 


± 1.4 


ACT-CL 


J2351. 7+0009 


1.5 


± 


0.2 


0.36 


± 


0.05 


1.5 


± 0.7 


3.2 


± 


1 





1.0 


± 1.4 


1.1 


± 


1.5 


3. 1 


± 1.6 


ACT-CL 


J0342. 7-0017 


3.4 


± 


0.3 


0.80 


± 


0.04 


3.2 


± 1.2 


3.1 


± 





.9 


1.1 


± 1.2 


5.0 


± 


1.5 


3.7 


± 1.4 


ACT-CL 


J0250.1+0008 


1.7 


± 


0.2 


0.41 


± 


0.05 


1.2 


± 0.5 


2.7 


± 





.8 


3.5 


± 1.1 


3.7 


± 


1.2 


3.0 


± 1.2 


ACT-CL 


J2152.9-0114 


1.9 


± 


0.2 


0.47 


± 


0.05 


1.5 


± 0.6 


3.0 


± 





.9 


3.9 


± 1.3 


1.2 


± 


1.1 


3. 1 


± 1.4 


ACT-CL 


J2130.1+0045 


1.9 


± 


0.2 


0.47 


± 


0.05 


1.6 


± 0.6 


3.2 


± 





9 


1.1 


± 1.2 


4.5 


± 


1.1 


3.7 


± 1.4 


ACT-CL 


J0018.2-0022 


1.8 


± 


0.2 


0.45 


± 


0.05 


1.5 


± 0.6 


3.1 


± 





.9 


1.0 


± 1.3 


1.1 


± 


1.1 


3.6 


± 1.4 


ACT-CL 


J0104.8+0002 


3.5 


± 


0.1 


0.81 


± 


0.06 


2.8 


± 1.3 


2.6 


± 





.9 


3.5 


± 1.2 


1.2 


± 


1.5 


3.0 


± 1.4 


ACT-CL 


J0017.6-0051 


1.6 


± 


0.5 


0.93 


± 


0.04 


5.0 


± 2.3 


2.9 


± 


1 





3.9 


± 1.4 


4.8 


± 


1.8 


3.4 


± 1.6 


ACT-CL 


J0230.9-0024 


2.5 


± 


0.3 


0.64 


± 


0.06 


1.9 


± 0.8 


2.8 


± 





.9 


3.7 


± 1.2 


4.3 


± 


1.1 


3.2 


± 1.4 


ACT-CL 


J0301. 1-0110 


2.0 


± 


0.3 


0.51 


± 


0.07 


1.1 


± 0.6 


2.2 


± 





,8 


2.8 


± 1.1 


3.2 


± 


1.3 


2.3 


± 1.2 


ACT-CL 


J0051. 1+0055 


1.7 


± 


0.2 


0.42 


± 


0.06 


0.9 


± 0.5 


2.2 


± 





.8 


2.8 


± 1.1 


3.0 


± 


1.3 


2.2 


± 1.1 


ACT-CL 


J0245.8-0042 


5.0 


± 


0.6 


0.95 


± 


0.04 


1.9 


± 2.4 


2.5 


± 





.9 


3.1 


± 1.3 


1.1 


± 


1.6 


2.8 


± 1.4 


ACT-CL 


J2051. 1+0056 


3.0 


± 


0.1 


0.73 


± 


0.07 


2.1 


± 1.1 


2.5 


± 





.9 


3.3 


± 1.3 


1.0 


± 


1.6 


2.8 


± 1.4 


ACT-CL 


J2135. 1-0102 


3.1 


± 


0.1 


0.75 


± 


0.06 


2.6 


± 1.2 


2.8 


± 


1 





3.7 


± 1.3 


4.5 


± 


1.6 


3.2 


± 1.5 


ACT-CL 


J0228.5+0030 


1.7 


± 


0.2 


0.42 


± 


0.06 


1.0 


± 0.5 


2.1 


± 





.8 


3.0 


± 1.1 


3.3 


± 


1.3 


2.5 


± 1.2 


ACT-CL 


J2135. 7+0009 


7.3 


± 


1.0 


1.01 


± 


0.01 


10.8 


± 6.3 


2.6 


± 


1 


.1 


3.6 


± 1.6 


1.3 


± 


1.9 


2.9 


± 1.8 


ACT-CL 


J2229.2-0004 


2.0 


± 


0.2 


0.49 


± 


0.06 


1.3 


± 0.7 


2.7 


± 


1 





3.1 


± 1.3 


3.8 


± 


1.5 


2.9 


± 1.4 


ACT-CL 


J0221. 5-0012 


1.6 


± 


0.2 


0.40 


± 


0.06 


0.5 


± 0.3 


1.4 


± 





.6 


1.8 


± 0.8 


2.0 


± 


0.9 


1.4 


± 0.7 


ACT-CL 


J2220.7-0042 


2.0 


± 


0.3 


0.50 


± 


0.07 


1.2 


± 0.7 


2.5 


± 


1 


.0 


3.1 


± 1.3 


3.5 


± 


1.5 


2.6 


± 1.4 


ACT-CL 


J2253.3-0031 


2.1 


± 


0.3 


0.54 


± 


0.06 


1.5 


± 0.7 


2.7 


± 





9 


3.5 


± 1.3 


4.0 


± 


1.5 


3.0 


± 1.4 


ACT-CL 


J2327.4-0204 


2.8 


± 


0.1 


0.67 


± 


0.01 


10.1 


± 1.0 


9.4 


± 


1 


,5 


12.5 


± 2.0 


14.3 


± 


2.2 


14.9 


± 3.0 



ACT: SZ Selected Galaxy Clusters 
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ID 0500 Q(m,z) Y 500 MUPP Mgg* M^ c 

(arcmin) (10~ 4 arcmin 2 ) (lO^/i^M©) (IO^/i^Mq) (10 14 /if 1 M G ) (lO 14 /^^©) 
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1 
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ACT-CL 


J2128.4+0135 


3.5 


± 





.1 


0.81 


± 





.02 


6.8 


± 


i.i 


5.7 


± 


i.i 


7.7 


± 


1 


6 


9.1 


± 


1 


.8 


8.3 


± 


2.1 


ACT-CL 


J2025. 2+0030 


3.6 


± 





.2 


0.83 


± 





.02 


5.5 


± 


1.3 


1.6 


± 


1.0 


6.3 


± 


1 


1 


7.7 


± 


1 


6 


6. 1 


± 


1.8 


ACT-CL 


J0026.2+0120 


2.2 


± 





.1 


0.56 


± 





.03 


2.8 


± 


0.6 


l.i 


± 


0.9 


5.8 


± 


1 


.2 


6.6 


± 


1 


3 


5.8 


± 


1.5 


ACT-CL 


J2307.6+0130 


3.4 


± 





.2 


0.80 


± 





.03 


1.5 


± 


1.2 


1.2 


± 


1.0 


5.7 


± 


1 


.3 


6.9 


± 


1 


6 


5.7 


± 


1.7 


ACT-CL 


J2156.1+0123 


1.9 


± 





3 


0.95 


± 





.02 


7.9 


± 
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1.1 


± 


1.0 


5.7 


± 


1 


1 


7.0 


± 


1 


7 


5.7 


± 


1.8 


ACT-CL 


J0301.6+0155 


6.7 


± 





1 
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± 





.01 


15.7 


± 


5.0 


1.7 


± 


1.2 


6.6 


± 


1 


.8 


8.1 


± 


2 


.1 


6.7 


± 


2.4 


ACT-CL 


J2051. 1+0215 


1.0 


± 





.3 


0.87 


± 





.03 


7.6 


± 


2.5 
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± 


1.4 


7.1 


± 


1 


.9 
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± 


2 


1 
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± 


2.6 
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14.9 


± 
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2.3 


ACT-CL 
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2.6 


± 





.2 
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.05 


2.1 


± 


0.8 


3.1 


± 


0.9 


4.0 


± 


1 


.2 


4.7 


± 


1 


1 


3.6 


± 


1.4 


ACT-CL 


J0219.9+0129 


3.0 


± 





.3 


0.73 


± 





.05 


2.3 


± 


0.9 


2.8 


± 


0.8 


3.7 


± 


1 


.1 


1.5 


± 


1 


1 


3.3 


± 


1.3 


ACT-CL 


J0240.0+0116 


2.1 


± 





.2 


0.53 


± 





.04 


1.8 


± 


0.6 


3.3 


± 


0.8 


4.3 


± 


1 


.1 


1.8 


± 


1 


.3 


3.9 


± 


1.3 


ACT-CL 


J0008.1+0201 


3.3 


± 





3 


0.78 


± 





.05 


4.0 


± 


1.6 


3.9 


± 


1.1 


5.2 


± 


1 


6 


6.2 


± 


2 





1.9 


± 


2.0 


ACT-CL 


J0139.3-0128 


1.6 


± 





2 


0.40 


± 





.06 


0.8 


± 


0.5 


2.1 


± 


0.9 


2.6 


± 


1 


.1 


2.8 


± 


1 


2 


2.1 


± 


1.1 



APPENDIX 



In this appendix we apply the methodology of Sections 3.1 and 3.2 to estimate masses and SZ quantities for the 
23 clusters from the original ACT Southern cluster sample, described in Marriage et al. (2011, hereafter Mil) and 
Menanteau et al. (2010). For this analysis, we use updated maps that include additional data acquired in this field 
during the 2009 and 2010 observing seasons. 

In Table 9 we summarize the basic properties of the Southern cluster sample, and include our measurement of the 
uncorrected central temperature decrement taken from maps filtered with the 6500 = 5f 9 UPP matched filter. While 
this filter incorporates the same cluster profile that was used in the analysis of the Equatorial sample, the filter itself 
is slightly different due to the different noise properties in the Southern maps. In Table 10 we present cluster mass 
and SZ quantity estimates. 

There is substantial overlap between the ACT Southern field and the cluster samples presented by the SPT collab- 
oration. The Mil sample includes six of the clusters analyzed in Williamson et al. (2011, hereafter Wll), and five of 
the clusters in Reichardt et al. (2012b, hereafter R12). In Wll a scaling relation based on the SPT signal to noise 
ratio and calibrated to the models of Shaw et al. (2010) is used to obtain SZ based mass estimates. In R12, masses 
are presented that are derived from a combination of SZ and X-ray measurements, though X-ray measurements are 
expected to dominate because of the smaller uncertainty in the X-ray scaling relation. (Four of the five ACT clusters 
appearing in R12 have X-ray mass measurements.) 

We find the SPT masses to be in good agreement with our masses based on the B12 model scaling relation parameters, 
with weighted mean mass ratio M 2 /M SPT = 0.99 ± 0.06. We note, however, a difference in the mean ratio for each 
of the two SPT catalogs. For the clusters in R12 the mass ratio is 0.88 ± 0.09, while for the higher mass, lower redshift 
sample of Wll the ratio is 1.10 ± 0.09. In Figure 21 we show the masses. We note that this difference in mass ratio is 
consistent with the observation in Wll that their average SZ determined masses are smaller than Yx based masses by 
a factor of 0.78 ±0.06. In contrast, the R12 masses for four of the five common clusters include input from observations 
of Yx, which is likely to dominate over the contribution from SZ data. 
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Figure 21. Comparison of ACT SZ based masses (from B12 scaling relation) to masses from the South Pole Telescope. Both are Msooc- 
ACT masses are computed using the B12 scaling relation and may be found in Table 10. SPT masses are taken from Williamson et al. 
(2011) (mass is computed from SZ signal based on scaling relation calibrated to Shaw et al. (2010) models; only statistical uncertainty is 
included in error bars), and Reichardt et al. (2012b) (mass is computed from combination of Yx an d SZ measurements). The weighted 
mean mass ratio for the 11 clusters is 0.99 ± 0.06, though there is evidence of a systematic difference between the mass calibration in 
the two SPT catalogs. Dotted line traces equality of SPT and ACT masses for the B12 scaling relation. Dashed lines trace approximate 
loci of agreement between SPT masses and ACT masses based on the UPP (upper line) and Nonthcrmal20 (lower line) scaling relations. 
Uncertainties in the ACT and SPT measurements for each cluster are likely to be partially correlated, since both use properties of cluster 
gas to infer a total mass. 



Table 9 

Table of uncorrected central Compton parameters for the clusters from 
Marriage ct al. (2011). Clusters used for cosmological and scaling relation 
constraints in Section 5 are marked with an asterisk. Redshifts are all 
spectroscopic, obtained from literature as described in Menanteau et al. 
(2010) or as presented in Sifon et al. (2012). We provide the S/N of detection 

as presented in Marriage et al. (2011), for reference. Uncorrected central 
Compton parameters yo are obtained from the application of a matched filter 

with 6500 = 5.'9 (see Section 3.1). Alternate IDs are provided for clusters 
detected before the initial release of the Southern sample, and for any cluster 
appearing in current SPT results. 





ACT ID 




R.A. 


Dec. 


Redshift 


S/N 






.W 






Alternate ID (ref.) 












(°) 


(°) 












(10- 


4 






ACT- 


-CL 


J0658- 


-5557 


104.6250 


-55.9511 





296 


11 


5 


2 


65 


± 





21 


lE0657-56/Bullet (1) 
































SPT-CL J0658-5556 


(2) 


ACT- 


-CL 


J0638- 


-5358 


99.6917 


-53.9792 





222 


10 





1 


77 


± 





15 


Abell S0592 (3) 
































SPT-CL J0638-5358 


(2) 
(2) 


ACT- 


-CL 


J0245- 


-5302 


41.3875 


-53.0344 





300 


9 


1 


1 


57 


± 





17 


Abell S0295 (3) 
SPT-CL J0245-5302 


ACT- 


-CL 


J0102- 


-4915* 


15.7208 


-49.2553 





870 


9 





3 


51 


± 





43 


SPT-CL J0102-4915 


(2) 


ACT- 


-CL 


J0438- 


-5419* 


69.5792 


-54.3181 





421 


8 





1 


62 


± 





13 


SPT-CL J0438-5419 


(2) 
(2) 


ACT- 


-CL 


J0645- 


-5413 


101.3750 


-54.2275 





167 


7 


1 


1 


19 


± 





17 


Abell 3404 (4) 
SPT-CL J0645-5413 


ACT- 


-CL 


J0546- 


-5345* 


86.6542 


-53.7589 


1 


066 


6 


5 





92 


± 





11 


SPT-CL J0546-5345 


(5) 


ACT- 


-CL 


J0235- 


-5121 


38.9667 


-51.3544 





278 


6 


2 





98 


± 





19 




ACT- 


-CL 


J0330- 


-5227* 


52.7250 


-52.4678 





442 


6 


1 


1 


23 


± 





18 


Abell 3128 NE (6) 




ACT- 


-CL 


J0616- 


-5227* 


94.1500 


-52.4597 





684 


5 


9 


1 


00 


± 





15 




ACT- 


-CL 


J0559- 


-5249* 


89.9292 


-52.8203 





609 


5 


1 





89 


± 





14 


SPT-CL J0559-5249 


(5) 


ACT- 


-CL 


J0215- 


-5212 


33.8250 


-52.2083 





480 


4 


9 





79 


± 





18 




ACT- 


-CL 


J0509- 


-5341* 


77.3375 


-53.7014 





461 


i 


8 





82 


± 





11 


SPT-CL J0509-5342 


(5) 


ACT- 


-CL 


J0641- 


-4949 


100.3958 


-49.8089 





146 


1 


7 





58 


± 





26 


Abell 3402 (4) 


ACT- 


-CL 


J0232- 


-5257 


38.1875 


-52.9522 





556 


1 


7 





61 


± 





17 




ACT- 


-CL 


J0516- 


-5430 


79.1250 


-54.5083 





294 


i 


6 





87 


± 





15 


Abell S0520 (3) 


































SPT-CL J0516-5430 


(5) 


ACT- 


-CL 


J0346- 


-5438 


56.7125 


-54.6483 





530 


i 


4 


1 


05 


± 





22 




ACT- 


-CL 


J0217- 


-5245 


34.2958 


-52.7556 





343 


i 


1 





79 


± 





18 


RXC J0217.2-5244 (7) 


ACT- 


-CL 


J0145- 


-5301 


26.2458 


-53.0169 





118 


1 








86 


± 





19 


Abell 2941 (4) 




ACT- 


-CL 


J0237- 


-4939 


39.2625 


-49.6575 





334 


3 


9 





93 


± 





26 





ACT: SZ Selected Galaxy Clusters 

Table 9 — Continued 
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ACT ID 




R.A. 


Dec. 


Redshift 


S/N 


yo 


Alternate ID (ref.) 








(°) 


(°) 




(io- 4 ) 


ACT- 


-CL J0304~ 


-4921 


46.0625 


-49.3617 


0.392 


3.9 


1.52 ± 0.32 




ACT- 


-CL J0707- 


-5522 


106.8042 


-55.3800 


0.296 


3.3 


0.54 ± 0.21 




ACT- 


-CL J0528- 


-5259 


82.0125 


-52.9981 


0.768 


3.1 


0.50 ± 0.13 


SPT-CL J0528-5300 (5) 



Note. — References: (1) Tucker et al. (1995); (2) Williamson et al. (2011); (3) Abell et al. (1989); (4) Abell (1958); (5) Reichardt ct al. 
(2012b); (6) Werner et al. (2007); (7) Bohringer et al. (2004). 

Table 10 

SZ-derived mass estimates for the ACT Southern cluster sample of 
Marriage et al. (2011). Columns are as described in Table 8. 



ID 05<x> Q(m,z) Y 5 oo Mgft, M^ c 

(arcmin) (10~ 4 arcmin 2 ) (10 14 h^M @ ) (10 14 /if M ) (l0 14 h^M Q ) (10 14 hf 1 M G ) 



ACT-CL 


J0658- 


-5557 


5.4 ± 


0.1 


0.97 


± 





01 


26.6 


± 


3.2 


10.3 


± 


1.9 


14.3 


± 


2 


.6 


17.6 


± 


3 


1 


18.2 


± 


1.3 


ACT-CL 


J0638- 


-5358 


6.1 ± 


0.2 


1.00 


± 





.01 


22.4 


± 


3.0 


7.5 


± 


1,1 


10.5 


± 


2 





13.0 


± 


2 


,1 


12.5 


± 


3.2 


ACT-CL 


J0245- 


-5302 


4.6 ± 


0.2 


0.92 


± 





01 


12.1 


± 


2.0 


6.7 


± 


1.3 


9.2 


± 


1 


,8 


11.3 


± 


2 


1 


10.4 


± 


2.6 


ACT-CL 


J0102- 


-4915 


2.5 ± 


0.1 


0.60 


± 





.02 


11.5 


± 


2.0 


10.5 


± 


1.8 


13.8 


± 


2 


,1 


15.0 


± 


2 


6 


16.0 


± 


3.5 


ACT-CL 


J0438- 


-5419 


3.6 ± 


0.1 


0.80 


± 





.01 


8.8 


± 


1.0 


7.0 


± 


1.2 


9.5 


± 


1 


.6 


11.5 


± 


1 


9 


10.9 


± 


2.3 


ACT-CL 


J0645- 


-5413 


6.9 ± 


0.3 


1.02 


± 





01 


18.2 


± 


4.3 


5.1 


± 


1.2 


7.3 


± 


1 


7 


8.9 


± 


2 


.0 


7.7 


± 


2,1 


ACT-CL 


J0546- 


-5345 


1.6 ± 


0.1 


0.36 


± 





.02 


2.0 


± 


0,1 


3.9 


± 


0.7 


5.1 


± 





.9 


5.2 


± 


1 





1.9 


± 


1.1 


ACT-CL 


J0235- 


-5121 


4.1 ± 


0.3 


0.87 


± 





.03 


5.9 


± 


2.0 


4.1 


± 


1.1 


5.6 


± 


1 


.5 


6.8 


± 


1 


.8 


5,1 


± 


1.9 


ACT-CL 


J0330- 


-5227 


3.2 ± 


0.1 


0.74 


± 





.02 


5.1 


± 


1.2 


5. 1 


± 


1.1 


7.3 


± 


1 


.5 


8.7 


± 


1 


7 


7.7 


± 


2.0 


ACT-CL 


J0616- 


-5227 


2.2 ± 


0.1 


0.53 


± 





.02 


3.0 


± 


0.6 


1.6 


± 


0.9 


6.1 


± 


1 


.2 


6.9 


± 


1 


3 


6.2 


± 


1.5 


ACT-CL 


J0559- 


-5249 


2.3 ± 


0.1 


0.56 


± 





.02 


2.8 


± 


0.6 


1.3 


± 


0.8 


5.7 


± 


1 


.1 


6.5 


± 


1 


3 


5.6 


± 


1,1 


ACT-CL 


J0215- 


-5212 


2.5 ± 


0.2 


0.62 


± 





.05 


2.1 


± 


1.0 


3.5 


± 


1.0 


1.5 


± 


1 


,1 


5.2 


± 


1 


7 


1.1 


± 


1.7 


ACT-CL 


J0509- 


-5341 


2.8 ± 


0.1 


0.66 


± 





.03 


3.1 


± 


0.7 


1.0 


± 


0.8 


5.3 


± 


1 


1 


6.3 


± 


1 


3 


5.2 


± 


1,1 


ACT-CL 


J0641- 


-4949 


4.8 ± 


1.0 


0.93 


± 





06 


2.9 


± 


3.2 


1,1 


± 


0.9 


1.8 


± 


1 


.2 


2.2 


± 


1 


,1 


1.3 


± 


1.1 


ACT-CL 


J0232- 


-5257 


2.0 ± 


0.3 


0.49 


± 





.07 


1.2 


± 


0.7 


2.5 


± 


1.0 


3.1 


± 


1 


3 


3.5 


± 


1 


.5 


2.6 


± 


1,1 


ACT-CL 


J0516- 


-5430 


3.9 ± 


0.2 


0.85 


± 





.03 


5.1 


± 


1.3 


1.0 


± 


0.9 


5.4 


± 


1 


.2 


6.7 


± 


1 


.5 


5.3 


± 


1.6 


ACT-CL 


J0346- 


-5438 


2.6 ± 


0.2 


0.62 


± 





.04 


3.3 


± 


1.1 


1.1 


± 


1.2 


5.7 


± 


1 


.6 


6.6 


± 


1 


.9 


5.5 


± 


2.0 


ACT-CL 


J0217- 


-5245 


3.2 ± 


0.3 


0.75 


± 





.05 


3.2 


± 


1.3 


3.1 


± 


1.0 


1.1 


± 


1 


.4 


5.3 


± 


1 


,8 


1.0 


± 


1.7 


ACT-CL 


J0145- 


-5301 


8.1 ± 


0.9 


1.02 


± 





.01 


16.9 


± 


8.2 


3.5 


± 


1.2 


4.8 


± 


1 


,8 


5.7 


± 


2 


.2 


1.2 


± 


2.2 


ACT-CL 


J0237- 


-4939 


3.2 ± 


0.5 


0.74 


± 





.09 


3.1 


± 


2.1 


3.1 


± 


1,1 


3.9 


± 


1 


.9 


4.6 


± 


2 


,1 


3.2 


± 


2.1 


ACT-CL 


J0304- 


-4921 


3.5 ± 


0.3 


0.79 


± 





.05 


6.8 


± 


2.6 


5.7 


± 


1.6 


7,1 


± 


2 


3 


8.7 


± 


2 


,8 


7.3 


± 


3.1 


ACT-CL 


J0707- 


-5522 


2.8 ± 


0.5 


0.69 


± 





.09 


1.1 


± 


1.1 


1.7 


± 


0.9 


2.2 


± 


1 


.2 


2.6 


± 


1 


4 


1.7 


± 


1.1 


ACT-CL 


J0528- 


-5259 


1.6 ± 


0.2 


0.37 


± 





.05 


0.9 


± 


0,1 


2.3 


± 


0.8 


2.9 


± 


1 





3.1 


± 


1 


1 


2,1 


± 


1.1 



